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Abstract 

We study electron-positron annihilations into six jets at the parton level in 
perturbative Quantum Chromo-Dynamics (QCD), via the elementary processes 
e + e _ —*• qqgggg, e + e~ —*• qqq'q'gg and e + e _ —> qqq'q'q"q" , for massive quarks 
q, q' and q" and massless gluons g. Several numerical results of phenomenological 
relevance are given, at three different collider energies and for a representative 
selection of jet clustering algorithms. We also present helicity amplitudes and 
colour factors needed for the tree-level calculation. 
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1 Introduction and motivation 


1.1 Multi-jet rates 

As the centre-of-mass (CM) energy of modern accelerators grows, the number of final state 
partons that can be accommodated in the available phase space becomes larger. Many of 
these will be quarks and gluons, as QCD forces are the strongest among the fundamental 
interactions. In electron-positron annihilations, coloured partons can be produced via an 
initial splitting of a 7 or Z current into quark-antiquark pairs, followed by no hard radiation 
at all or by (multiple) gluon bremsstrahlung, with the vector particles eventually yielding 
further quark-antiquark pairs or other gluons, as schematically recalled in Fig. []]. 



etc. 


Figure 1: Schematic description of multi-parton production in e + e annihilations. 

Clearly, as this so-called parton shower or cascade evolves in time, the amount of energy 
exchanged among the various partons diminishes, so that the process will eventually enter a 
non-perturbative regime. At such low energy scales (roughly, of the order of 1 GeV), one is 
unable to describe quantitatively the above process in terms of a series expansion in powers 
of the strong coupling constant a s using Feynman diagram techniques, since a s becomes 
significantly large and the series cannot be resummed reliably. Note that 1 GeV is also 
the scale at which the partons are conventionally assumed to hadronise, that is, to convert 
themselves into the observed hadrons through the formation of colourless bound states (see 
Ref. |I|] for an up-to-date review of hadronisation mechanisms), which in turn make up the 
jets observed in the experimental apparatus. 

This fact clearly illustrates an intrinsic limitation of the perturbative approach. A fur¬ 
ther, more practical, restraint that one has to face is due to our scarce ability for calculating 
multi-particle amplitudes using standard perturbative techniques. I 11 fact, with the number 


of final states particles increasing, the complexity of the expressions of the fundamental scat¬ 
tering amplitudes and the large number of terms involved (arising from the combinatorics 
of the diagrams) that one has to master become in the end an obstacle difficult to overcome, 
in spite of the help of modern computers. For example, if one limits oneself in Fig. [I] to the 
case of gluon radiation only (that is to qq + ng final states), one has the following impressive 
escalation of numbers of Feynman diagrams involved (considering the 7 and Z propagators 
together): 1 for n = 0, 2 for n — 1, 8 for n = 2, 50 for n = 3, 428 for n = 4, 4670 for n = 5, 
etc]]. 

An even higher degree of complications is encountered if one goes beyond the leading 
order (LO) and tries to incorporate loop corrections in the calculation. In fact, it should be 
recalled that the actual number of jets produced in electron-positron annihilations is defined 
according to a phenomenological jet clustering algorithm, in which jets are constructed out of 
primary objects (hereafter, labelled by i,j etc.), the latter being hadrons or calorimeter cells 
in the real experimental case and partons in the theoretical calculation. The term exclusive 
means that each primary object is assigned to a unique jet and each final state has a unique 
jet multiplicity, for a given value of the jet resolution parameter y cut . Algorithms of this 
type are generically based on a binary clustering, wherein the number of objects is reduced 
one at a time by combining the two nearest (in some sense) ones. The joining procedure 
is stopped by testing against some criterion (the parameter y cut ), and the final clusters are 


called jets (we will treat some of these algorithms in greater detail later on, in Subsect. | 2 . 2 | ). 

From the point of view of perturbative calculations, it is evident that the ?r-jet rate 
will then be constituted not only of the fraction of n-parton events in which all partons 
are resolved, but also by the (n + m)-parton contributions (with m > 1 ) in which m 
partons remain unresolved. Now, in the (n + m)-parton —> n-jet transitions, two partons 
can well be combined in singular configurations. This can happen when one has either two 
collinear partons within or one soft parton outside the ‘jet cones’. Both these contributions 
are (in general, positively) divergent. Thanks to the Bloch-Nordsieck [(]] and Kinoshita-Lee- 
Nauenberg § theorems (see also Ref. 0), these collinear and soft singularities are cancelled 
at the same order in a s by the divergent contributions (generally, negative) provided by the 
virtual loop graphs. If one recombination has occurred then one-loop diagrams must be 
computed, when two mergings have taken place two-loop graphs are needed, and so on. As 
a matter of fact, the computation of loop diagrams is even more demanding than that of 
tree-level ones (both using trace [7J and helicity techniques). Not surprisingly then, the 
number of exact calculations of e + e~ —> n-parton matrix elements (ME) available in the 

2 Under these circumstances, the use of spinor techniques based on the helicity decomposition of the 
Feynman amplitudes becomes mandatory, in order to simplify the calculation to a certain extent (for a 
good introduction to helicity methods and complete references, see 0 and S). 




literature is rather limited. 


As for the tree-level, till a few years ago only up to hve-parton final states were computed, 
via the elementary scatterings: e + e“ —> qqggg and e + e _ —» qqq'q'g. This calculation was 
tackled in Ref. 0. The case of four partons via e + e _ —> qqgg and e + e“ —> qqq'q' has been 
known since long time [|1(J. The case of three partons (i.e., e + e~ —> qqg) was very trivial 
||T[ and the lowest-order case of two (i.e., e + e _ —> qq) is in fact a textbook example [|TT . 
Only recently, the task of computing six-parton cross sections (via e + e _ —> qqgggg, e + e~ —> 
qqq'q'gg and e + e~ —> qqq'q 1 q"q") has been accomplished, in Ref. |H|, where numerical 
results for the massless case were presented. It is the main purpose of the present paper 
to present new phenomenological analyses of six-jet production at lepton colliders as well 
as to give the analytical expressions corresponding to those six-parton processes, further 
extending the calculations of Ref. [jl2) to the case of quarks with generic mass. Notice 
that the retention of finite values of m qtq > tq n in the amplitudes turns out to be a serious 
complication in many instances, because of the unavoidable proliferation of terms and/or 
the slowness of execution of the numerical routines that it generates. As a matter of fact, 
depending on the CM energy, the actual flavour of the quark being produced and also the 
region in phase space under study, in many cases the effect of the masses can safely be 
ignored^. 

As for loop and radiative diagrams, the highest-order corrections calculated to date in 
e + e _ —» rz.-jet annihilations are the next-to-leading order (NLO) ones to the four-jet rate, 
that is, terms proportional to [[H], 0, in the massless approximation. The ones 
date back to Ref. m for the massless case while the massive one has been completed very 
recently [[Rj. In fact, considerable simplifications can be attained also at loop level when 
the quark mass is neglected. Two- and higher-order-loop corrections are not yet known. 

But let us review the theoretical progress in evaluating e + e” —> n-parton rates and its 
interplay with the experiment. This will also eventually help us to motivate the need of the 
availability of the calculations carried out in Ref. |jl2| and in this paper. 

The history of calculations of multi-jet production in e + e~ scatterings |lf|-[|28| began 
at the times when jets were first observed at PEP and PETRA. At those energies the Z 
contribution was small and could safely be ignored ||TD| , [P^H~|I5|, (see also JH], j2H|). Mass 
effects were already taken into account in some instances. However, apart from the simplest 
cases [00, only the formulae for massless fermions were compact enough to be published 
mmm- With the advent of LEP and the advances in spinor techniques these results were 
improved in various ways. In particular, in |25|] the complete triply differential cross section 


3 So was the case before the LEP era for the more complicated contribution of the intermediate Z (and 
its interference with the photon) in the primary e + e _ annihilation, provided the energy at which the latter 
takes place is well below the mass of the former (at PEP and PETRA, for example). 





















































for e + e~ —» qqg for massive quarks with both intermediate photon and Z was presented. In 
Ref. [^] the decay Z —> 4/ (with f = q or £: i.e., quark or lepton) was studied and from the 
expression given there it is possible to reconstruct the full amplitude for e + e _ —> qqq'q '. Very 
simple formulae for e + e~ —> qqgg , qqq'q 1 were given in Ref. |27]] for massless quarks. The 
tree-level amplitudes for three-, four- and five-parton events including all mass effects and 
both the 7 and Z contributions were given in analytic form in Refs. |29], |3(J. Finally, one- 
loop three-parton and tree-level four-parton terms have been combined and incorporated in 
Monte Carlo simulation programs |R], |32], [k§ ; as has been done for the one-loop four-parton 
and tree-level five-parton amplitudes [[T| 


36| , for massless quarks in both cases, whereas 
the massive contributions to the aj, three-jet rates with massive partons will soon be. 

A further point should be noted. On the one hand, at small values of the resolution 
parameter y cut , the fixed-order perturbative calculations of n-jet rates become unreliable 
because of the presence of large logarithms of the form ln(y cut ) (in other words, because of 
the growing of the (n + m)-jet contributions, with m > 1). On the other hand, the exact 
higher order calculations can be technically unavailable. Under these circumstances, the 
resummation to all a s orders in perturbation theory of terms of the form 0(a™ ln m (l/y)), 
e.g., with m = l and m — l — 1 (to leading and next-to-leading accuracy, with n and m 
depending on the jet observable), can improve substantially the reliability of the theoretical 

Ell 


predictions, as was shown to be the case for the two- and tree-jet rates 


Recall, 


however, that such a resummation is only possible within the framework of certain jet 
clustering schemes, such as the Durham and Cambridge ones (see Subsect. |2]2|), for which 
the jet fractions show the usual Sudakov exponentiation of multiple soft-gluon emission. 

The availability of three-, four- and five-parton matrix elements (tree- and loop-level and 
supplemented with the mentioned logarithms resummed) has been of crucial importance in 
studying the underlying structure of QCD. Throughout the years and at different colliders 
(PEP, PETRA, LEP, SLC), several tests have been performed. The strong coupling constant 


a s has been determined from jet rates and from shape variables [42| and both the flavour 
independence and the running with y/s have been verified (see, e.g., Ref. [[T|). Three- 
[[44l and four-jet (45| differential distributions have been long studied and their behaviour 
agrees with QCD predictions calculated to second and third order in ct s , respectively. The 
colour factors, which determine the gauge group responsible for strong interactions, have 
been measured |R| pT7| . Models alternative to QCD have been ruled out and the coupling 
of the triple gluon vertex has been verified to be consistent with the QCD theory |I3|] . The 
latest QCD tests are currently being carried out at LEP2. In addition, at this energy, one 
of the key issues is to quantify to a high degree of accuracy the amount of background 
in multi-jet events from QCD to the W + W~ hadronic decays, as the corresponding signal 
(though biased by systematic uncertainties due to relatively unknown Bose-Einstein El 



































and Colour Interconnection [j50| effects) still represents to date the experimentally preferred 
decay channel, because of its somewhat higher statistics and because kinematic constraints 
can tighten the precision of M w ± measurements [[51], |52|| . 

In the not too distant future (presumably not long after the year 2005), a new electron- 
positron machine will be operative, the Next Linear Collider (NLC), which is expected to 
run initially at an energy of 500 GeV, later increasing up to 2 TeV or so. 

1.2 Six-jet rates 

The phenomenology of six-jet events produced in electron-positron annihilations (shown 
schematically in Fig. j2j) is not as all well known as that of smaller multiplicity jet rates, 
from both the experimental and theoretical point of view. This should not sound surprising, 
given the complications arising from, on the one hand, the poor event rate and large number 
of tracks, and, on the other hand, the mentioned complexity of the perturbative QCD 
calculations. 
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Figure 2: Schematic Feynman diagrams contributing in lowest order to e + e _ —> 6 partons. 
The jagged and dashed lines and the labels x, y refer to either a quark or a gluon, see (|])-(|3|). 

In our opinion, the time has now come to approach the calculation of e + e _ —> 6 jet 
events. We are driven to this belief from the following considerations. LEP1 was the era 
of the Z-peak and of its two-jet (and several higher-order) decays. LEP2 is now the age 
of the iy + IE - -resonance and of its four-jet (and some higher-order) decays. As we will 
move into the NLC epoch, we will step into a long series of resonant processes, typically 
ending up with six-jet signatures. One can mention top physics for a start. As this is one 
of the main goals of the NLC, we should expect a lot of experimental studies concerned 
with tt —> bbW + W~ —> 6 jet-events, since the W ±, s mainly decay hadronically into two 
jets. Then one should not forget the new generation of gauge boson resonances, such as 
W + W~Z and ZZZ. and their favourite six-jet decays. One could also add highly-virtual 



photonic processes, like ^W + W~ , 7 ZZ, 77 Z and 777 as well as those involving the Higgs 
particle, via ZH —> ZkH + IR _ , Zf/ —»■ ZZZ and Zf/ — > ZHH, provided this is light enough 
to be produced at those energy regimes. 

A collection of dedicated articles illustrating the phenomenology of top-antitop and 
three-boson production at NLC can be found in Ref. |53|. It is beyond the scope of this 
paper to recall the salient features of those studies and the interplay of the mentioned 
resonant events with the six-jet background produced via QCD. Instead, what we would 
like to emphasise here, is that one should be ready with all the appropriate phenomenological 
instruments to challenge the new kind of multi-jet experimental studies that will have to be 
carried out at the new generation of e + e~ machines. 

As for the theoretical progress in this respect, studies of e + e“ —> 6 -fermion electroweak 
processes are under way. A brief account of approaches and methods can be found in 
Ref. |>|]. These kinds of reactions have already been calculated and analysed for the case 
involving up to four quarks in the final state ||55| (and other similar studies are in preparation 
fjEf), while in Ref. |57J the concern was for Higgs processes with two quarks produced. One 
should probably expect the upgrade of the codes used for those studies to the case of 
six-quark production via electroweak interactions quite soon. However, a large fraction of 
the six-jet cross section comes from QCD interactions involving gluon propagators, gluon 
emissions and quark-gluon couplings. The case of six-jet production from W + W~ decays 
was considered in Ref. [ |58| . 

In the present paper, the main focus will be on the 0(af) QCD processes at tree-level: 


e + e 

QQ9999, 

( 1 ) 

e + e _ ■ 

-*• QQQ'Q'99 

( 2 ) 

• + e~ - 

- / -/ // -// 

+ QQQQQ Q , 

(3) 


where q,q' and q" represent any possible flavours of quarks (massless and/or massive) and 
g is the QCD gauge boson. 

As a by-product of our analysis, we will study the case of electromagnetic (EM) radiation. 
Since the QED emission of photons off quark lines is dynamically similar to that of gluons 
(apart from the different strength of the coupling and the colour factor associated with it), 
we will modify our MEs with external gluons (|T])~(@) appropriately, so as to allow for the 
case of multiple hard photon emission from the final state (Final State Radiation, FSR), in 
e + e _ —*• qqggg'f, e + e~ —> qqgg'y'y, e + e~ —> <717(7777, e + e~ —> <7(77777 (via two-quark) and 
e + e _ —> qqq'q’gj, e + e~ —> qqq'q '77 (via four-quark) events. I 11 order to calculate the pho¬ 
toproduction cross sections, we further supplement the FSR diagrams with those in which 
photon emission takes place from the initial state (Initial State Radiation, ISR) along with 







those in which photons are radiated from both the incoming and outgoing fermions. Notice 
that hard photons produced by electrons and positrons can also be included by means of a 
convolution of the non-radiative diagrams with the so-called Electron Structure Functions 
(ESF) j59|. Such implementation allows one to include the exact photon corrections to the 
e + e~ annihilation subprocess resummed up to the third order |Kj in the EM coupling con¬ 
stant a em in the infrared limit (in particular, they embody both soft/collinear and virtual 
photon emission). We have not exploited here this approach, as we have used our tree-level 
diagrams. In doing so, however, additional cuts are necessary in order to screen our results 
from the collinear singularities due to the infrared emission from the initial state. We have 
required p? > 5 GeV and cos 0 7 1 < 0.85. Such cuts are similar to the photon selection 
criteria adopted by the LEP Collaborations. 

We will apply our results to the case of LEP1, LEP2 and, particularly, NLC energies. A 
short description of the computational techniques adopted will be given in Sect. |^. Results 
are in Sect. [3]. A brief summary is given in Sect. f|. The explicit helicity amplitude formulae 
are in Appendix I, whereas the colour factors are presented in Appendix II. 


2 Matrix Elements 


In order to master the large number of Feynman diagrams (of the order of several hundred) 
entering in processes (|I])-(|3|) at tree-level, we have used spinor methods [|3|. In fact, if one 
adopts the technique of deriving the expression of the squared amplitude \A\ 2 = \ Y.k a k\ by 
using the textbook method of the traces (a*, being the amplitude associated to diagram k ), 
one results in a prohibitive number of terms to be evaluated. Typically, if n is the number 
of diagrams, then n A+ l ) traces have to be computed. Recalling the numbers of Sect. jl]. 
the task is clearly prohibitive for six-parton events. The calculation becomes simpler if one 
uses so-called ‘helicity amplitude’ (or ‘spinor’) techniques. In such an approach, given a 
set of helicities assigned to the external particles, one computes the contribution of every 
single diagram k as a complex number a*,, sums over all /c’s and takes the modulus squared. 
This way, in the above example, one only needs to compute n (complex) terms. To obtain 
unpolarised cross sections one simply sums the modulus squared for the various helicity 
combinations. Thus, a further advantage is intrinsic to the spinor method, that one can 
naturally compute also polarised MEs (e.g., for the SLC and NLC), without the need of 
inserting into the expression of the amplitudes the helicity projection operators |JI . 









2.1 The spinor methods 


The use of helicity amplitude techniques in high energy physics dates back to Refs. |[(T2|. F3 


Since then, many different approaches have been developed J3|-[|75| and a review of the 
vast literature on this subject is beyond the intentions of this paper. For this, we remind 
the reader of Ref. [J2|. Of those available in literature, we employed here two methods. We 
have done so in order to be able to check adequately the correctness of our calculations, 
as usual tests of gauge invariance might fail to reveal possible bugs in the numerical codes. 
Note that the two methods can be used both with massless and massive particles. 

The first one is based on the formalism of Ref. |72| and exploits the HELAS subroutines 
It is the same already employed to produce the numerical results presented in [|12]. 


76 


Here, one adopts a specific representation for the Dirac matrices. It is then possible to 
obtain explicit expressions for the spinors, for the matrices at the vertices and for the fermion 
propagators, and thus the components of the polarisation vectors of the gauge bosons can 
be written directly in terms of their four-momenta. The complex number corresponding to a 
given diagram is obtained by multiplying these matrices. Therefore, one most conveniently 
programs the amplitudes directly as deduced from the Feynman rules and no manipulations 
on these need be performed. Masses in this case turn out to be an entry in a propagator 
matrix. It would be possible to consider some algebraic simplifications of the formulae, 
e.g., by separating the two independent chirality components of the massless spinors and 
cancelling by hand the terms proportional to the quark masses. However, the method of 
Ref. |72j is particularly well suited for numerical calculations and such algebraic treatment 
would not significantly improve the speed of the programs. Therefore, we have not exploited 
here any chiral decomposition in case of massless spinors. Under these circumstances, the 
analytical expressions in the formalism of Ref. |72| would simply be a rewriting of the 
Feynman rules with the mentioned matrices expressed explicitly in terms of the external 
momenta and hclicities, an exercise that we do not perform here. 

A second approach was based on the method developed in Refs. [7]. [73|, |74], [77|| . Since 
we will express the helicity amplitudes in this formalism, we will devote same space to its 
technicalities. (This is also helpful in order to introduce a notation that will be used in the 
Appendices and to render the formulae appearing there more transparent to the reader.) 

In this method, all spinors for any physical momentum are defined in terms of a basic 
spinor of an auxiliary light-like momentum. By decomposing the internal momenta in terms 
of the external ones, using Dirac algebra and rewriting the polarisation of external vectors 
by means of a spinor current, all amplitudes can eventually be reduced to an algebraic 
combination of spinors products u(pi)u(pj) (with i,j ,... labelling the external particles). 






































A 1 A 3 

X(p 1 ,Xi;p 2 ;p 3: X 3 ]c R ,c L ) 

++ 

(p1V2 + P 2 Vi)(crP 2 V 3 + c L p 3 p 2 ) + c R S(+,p 1 ,p 2 )S(-,p2,p 3 ) 

+- 

C L (P 1V2 + P 2 Vl)S(+,P 2 ,P 3 ) + C L (c L p 2 V 3 + C R p 3 ?) 2 )S( +, p u p 2 ) 


Table 1: The A" functions for the two independent helicity combinations in terms of the 
functions S, ij and p defined in the text. The remaining X functions can be obtained by 
flipping the sign of the hclicities and exchanging + with — in the S functions and R with 
L in the chiral coefficients. 


A1A2A3A4 

^{Pl , Ai, P2') A 2, P3i A3, P4i A 4 , C R , C L , Cfo C L ) 

+ + ++ 

-2[S(+,p 3 ,p 1 )S(-,p4,p 2 )c , R c R - piP 2 imAc' R c L - rjirj2P3P^ L c R ] 

+ + +- 

-2-// 2 C i? [^( + ,p 4 ,Pl)/i3Ci - S(+,P3,Pi)P4c' r ] 

+ + -+ 

-27] 1 C R [S(-,P2,P 3 )P4C' L - S(-,P2,P 4 )P3C' R ] 

+ - ++ 

-27]4c' R [S( + ,P3,pi)p 2 C R - S(+,P3,p 2 )piC L \ 

+ + — 

—2{S(+, Pi, p 4 )S( —, p 2 , p 3 )c' L c R - P 1 P 2 V 3 Inc' L c L - rjirj 2 p 3 P4^ R c R ] 

+ - +- 

0 

+ - -+ 

—2 \pip 4 r] 2 il 3 c' L c L + P 2 P 3 V 1 V 4 C R C R - P 2 P 4 V 1 V 3 c' L CR - P 1 P 3 V 2 V 4 d R c L \ 

+ - - . 

-2t]3c'l[S(+ 1 P2iP4)PiC L - S(+,p 1 ,p 4 )p 2 c R } 


Table 2: The Z functions for all independent helicity combinations in terms of the functions 
S, 7] and fi defined in the text. The remaining Z functions can be obtained by flipping the 
sign of the hclicities and exchanging + with — in the S functions and R with L in the chiral 
coefficients. 


(i) Spinors. External fermionsQ of mass m and momentum are described by spinors 
corresponding to states of definite helicity A, u(p, A)Q verifying the Dirac equations 

y(u(p, A) = ±mu(p, A), u(p, AW = ±mu(p,\), (4) 

and the spin sum relation 

U (P> X MP, A) = / ± m, (5) 

A=± 

where the sign +(—) refers (here and in the following) to a particle(antiparticle). One can 

4 Unless stated otherwise, we shall use the term “fermion” and the symbol “u” for both particles and 
antiparticles. 

5 Here, p( A) represents a generic (anti)spinor four-momentum(helicity). 



















choose two arbitrary vectors ko and k\ such that 


k 0 ■ k 0 = 0, ki ■ k\ — —1, ko ■ ki — 0, (6) 

and express the spinors u(p, A) in terms of chiral ones w(k 0 , A) as 



u(p, A) = w(p, A) + nw(k 0 , -A), 

(7) 

where 

w(p, A) =/w(k Q ,~ A)/?7, 

( 8 ) 

and 

m j— -— 

/i = ±—, p = J2\p-k 0 \. 

r/ v 

(9) 

The spinors w(k 0 , A) satisfy 

w(k Q , X)w(k 0l X) = 2 

(10) 

and therefore 

w(k Q , X)w(k Q , X) = 

(11) 


A=± 

The phase between chiral states is fixed by 


w(k Q , A) = \lf w{k Q , -A). 


( 12 ) 


The freedom in choosing k 0 and k\ provides a powerful tool for checking the correctness of 
any calculation. A convenient, though not unique choice, is the following: k 0 = (1, 0, 0, —1) 
and k\ = (0,1, 0, 0) [[74). In such a case the massless spinors in the two methods and m 


coincide exactly, so that it is possible to compare in greater detail the two corresponding 
numerical codes. In particular, the results obtained with the two formalisms must agree for 
every single diagram and every polarisation of external particles. 


(ii) Polarisation vectors for massless gauge boson^_ ]. External spin 1 massless gauge 
bosons of momentum p M are described by polarisation vectors corresponding to states of 
definite helicity A, £^(p, A), fulfilling the identities 


£(P, X) - p = 0, 

£(P, X) ■ e(p, A) = 0, 

(13) 

^(P,~X) =£ /1 *(p,X), 

£{P, A) • e(p,-A) = -1, 

(14) 


6 Polarisation vectors for massive gauge bosons can be introduced on a similar footing, however, since we 
do not need them for our calculation, we refer the reader to the mentioned bibliography for this particular 


case. 







and the spin sum relation (e.g., in the axial gauge) 


E ^(p, A)e-(p, A) = -<r + q ' Y + qUp \ 
Ati P-9 


( 15 ) 


where (f is any four-vector not proportional to p+ Any object £^{p, A) obeying the relations 


(|T4|)—(|T5D makes an acceptable choice for the polarisation vectors. For instance 

P, A) = N[u{p, \)Yu(q, A)], (16) 

N being the normalisation factor 

lV=[4(g-p)]- 1/2 . (17) 


The existing freedom in choosing corresponds to fixing the gauge. The final results do 
not depend on the choice of q A 


(iii) The S, X and Z functions. Using the above definitions one can compute the 
functions 

S(X,Pi,P2) = [u(pi, A)u(p2, -a)], (18) 

X(p 1 ,\ 1 -,p 2 -,p 3 ,\ 3 ;c R ,c L ) = [u(pi, Ai)|4 Tu(p 3 , A 3 )], (19) 

and 

^ ijp 1 : A 4 , P2j A 2 , PS: A 3 , P 4 , A 4 , Cr, Cl, Cr, (‘j ') 


where 


and 


with 


[u(Pu Ai)PXp 2 , A 2 )][u(p 3 , A 3 )r^«(p 4 , A 4 )], 


p(') _ J!) p 1 J!) p 
1 — C R r R + C L I Jj , 




1+75 


Pr = 


75 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


the chiral projectors. 

By computing the resulting traces one easily finds (e 0123 = 1 is the Levi-Civita tensor) 

mm 


\ Q (pi ■ k o)(P2 ■ h) - (pi ■ h)(p2 ■ ko) + le^pXoKp&Z /q yj\ 

z{+,pi,p 2 ) — ^ -, (+4J 

V 1 V 2 

for the S functions and the expressions listed in Tabs. [L] and |2] for the X and Z functions, 
respectively. For the S functions, one has S(—,pi,p 2 ) = S(+,P 2 ,Pi)*, while the remaining 
X and Z functions can be obtained as described in the captions of Tabs. [T] [|. 








The colour factors have been calculated with two different methods, as illustrated in 
Ref. [|78[ and by means of the techniques reported in Ref. [[79]. Some possible expressions 
of both the amplitudes and the colour matrices will be given in the Appendices. The coded 


hclicity amplitudes have been checked for gauge invariance and integrated using VEGAS [80 


in both formalisms. The agreement at the (squared) amplitude level is always within 12 
significant digits in REAL*8 precision. 

Although the actual number of Feynman diagrams needed to compute processes (||])-(0) 
is very large, only a much smaller portion of different topologies[] describe these. In fact, in 
the case of two-quark-four-gluon diagrams the topologies are those shown in Fig. 0(a), eight 
in total. For four-quark-two-gluons, Fig. 0(b), one has ten and for six-quarks, Fig. 0(c), 
three topologies. The same implementation of each of the topologies pictured in Fig. 0 
can then be used several times by means of recursive permutations of the momenta of the 
external particles. Furthermore, each common ‘sub-topology’ to those in Fig. 0is saved and 
then reused in the same numerical evaluation. The actual number of Feynman diagrams 
corresponding to the topologies in Fig. 0, including all possible permutations of momenta 
and counting as either one or three each diagram involving a four-gluon vertex (which carries 
in fact three different colour structures), is given in Tab. 0. 


2.2 Jet clustering algorithms 

As jet clustering algorithms we have adopted the Jade (J) |0Tf , Durham (D) [00] and Cam¬ 
bridge (C) 00] scheme^. The first one uses the expression 

y j 2E i E j {l-cosO ij ) 

J s 

for the ‘jet measure’, whereas the other two are based on the same test variables 

d c 2 min (Ef, Ej)(l - cos %) 

Vij = ---, (26) 

s 

but differ in the clustering procedure. In eqs. (]J5|)- (|26|) , E % and Ej represent the energies of 
any pair of partons i and j, whereas 9 l3 is their relative angle. A six-jet sample is selected 
by imposing the constraint y l3 > y on all n possible parton combinations ij (with n =15 

' By ‘topology’ we mean a collection of Feynman diagrams which differ only in the permutations of gauge 
vectors along the fermion lines (neglecting their colour structure). 

8 In this Section as well as in the following ones, in order to simplify the notation, we shall use y to 
represent y cut , the jet resolution parameter discussed in Sect. [T[ In addition, we acknowledge here our abuse 
in referring to the various jet ‘finders’ both as algorithms and as schemes, since the last term was originally 
intended to identify the composition law of four-momenta when pairing two clusters. This is in fact a well 
admitted habit which we believe will not generate confusion in our discussion. 














e+e -> qqgggg 

al 

a 2 

a3 

a4 

a5 

a 6 

a7 

a 8 



120 

144 

24 [72] 

72 

36 

8 [24] 

12 

12[36] 



Total 428 [516] 

e+e -> qqq'q'gg 

(q q'] 


bl 

b 2 

b3 

b4 

b5 

b 6 

b7 

b 8 

b9 

blO 

12 

48 

48 

24 

24 

8 

4[12] 

8 

4 

16 

Total 196 [204] 

e+e^ —> qqq'q'q”q" (q ^ q' ^ q") 

cl 

c 2 

c3 








24 

18 

6 








Total 48 


Table 3: Number of Feynman diagrams corresponding to the topologies in Fig. |2j. If two 
quark flavours are identical, numbers are multiplied by two, if three are, by six. In square 
brackets are the same numbers if the four-gluon diagrams are counted thrice, according to 
their three colour structures. 


for the Jade and Durham and 5 < n < 15 for the Cambridge algorithm). The n = 15 pairs 
of the J and D schemes are simply obtained by pairing all partons in the final state in all 
possible different combinations ij . The variable number 5 < n < 15 for the C scheme is a 


consequence of the so-called ‘soft freezing’ step, which is described in [|3|], to which we refer 
the reader for further details. 


2.3 Numerical parameters 

We have adopted Mz = 91.17 GeV, Y z = 2.516 GeV, sin 2 6w = 0.23, a em = 1/128 and, 
unless otherwise stated, the two-loop expression for a s , with Aqcd 4 = 0.298 GeV. We have 
kept all quarks massless as a default, to speed up the numerical evaluation. Though, in 
several instances we have retained a finite value for the bottom mass, i.e., = 4.25 GeV. 

Electron and positron have mass zero too. When e+e - —> tt —> bbW + W~ —> 6-jet events 
have been calculated (at tree-level), the following values of t- and VF+-parameters were 
adopted: m t = 175 GeV, r t = 1.54 GeV, M w = 80.23 GeV and Y w = 2.2 GeV. As centre- 
of-mass (CM) energies representative of LEP1, LEP2 and NLC, we have used the values 
a fs = M z , 180 GeV and 500 GeV, respectively. 

The vertex couplings (or ‘chiral coefficients’) Cr and Cl of eq. (|22|) between 
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z 

W 

9 

Cr 

Q f 

9r/ s wCw 

0 

1 

c l 

Q f 

9l/ s wCw 

1/\/2% 

1 


Table 4: The couplings Cr and Cl of eq. (|2j) for u- and d-type (anti)quarks and elec¬ 
trons/positrons to the gauge bosons 7 , Z, W and g. One has (adopting the notations 
sw = sin 0\v and cw = cos 6w) 9r = ~Q *and y{ = T/ — (with q — u, d), where 

(Q u , T 3 “) = (+|, +|), ( Q d , T 3 d ) = (-|, —|) and (Q e , Tf) = (-1, —|) are the fermion charges 
and isospins. 


quarks/leptons and gauge bosons entering in the expressions of the S', Y and Z functions 
can be found in Tab. £| (we have included also the W ± qq' vertex, for completeness, though 
it will not appear in our formulae). 


3 Results 

We define the ^/-dependent six-jet fraction by means of the relation 

* t N (d) Mu) fr>n\ 

Mv) = ^ ^ 7 x = ——, (27) 

zZm&myy) 

where cr 6 (y) is the actual six-parton cross section and a t identifies the total hadronic rate 
(T t = Co(1 + OLsj-n + ...), <T 0 being the Born cross section. I 11 perturbative QCD one can 
rewrite eq. (|27|) in terms of a series in ck s , beginning with its fourth power, as 


Mv) - (2^) My) + •••> 


(28) 


where G(y) is the lowest-order (or leading-order, LO) ‘coefficient function’ of the six-jet rate. 
We show this quantity in Fig. [|. for the case of the Durham, Cambridge and Jade schemes, 
at LEP10. Our rates for the Jade and Durham schemes nicely match the ‘leading-order’ 
predictions (both in a s and log(l/y)) given in Ref. |84|] in the region of y ~ 1 x 10 ~ 4 for 
the Durham scheme and when y & 5 x ICR 4 for the Jade one. Here and in the following, 
though, we will take as a lower limit of our y-range the value of 0 . 001 , where fixed-order 
perturbation theory should be much more reliable. 

A salient feature of Fig. [f]is that the Jade rate is much above the Durham and Cambridge 
ones. At y = 0.001 the ratio between the former and the latter two is about 8.0. (Notice 

9 Here and in the following, unless otherwise stated, the summations over the three reactions (0) 1) and 
over all possible massless combinations of quark flavours in each of these have been performed. 
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Figure 4: The parton level G(y) function entering in the six-jet fraction at LO in the D 
(continuous line), C (dashed line) and J (dotted line) schemes, at LEP1. 


that the corresponding number for the three-, four- and five-jet rates is 2.3, 3.6 and 5.7, 
respectively.) This effect is due to the fact that y'- 3 is always greater than or equal to y ®’ c , 
as can be easily understood from eqs. fl25l) -([26l) . Therefore, for a specified y and a given 
configuration of energies and angles E i: Ej , cos 6 tJ , more events will pass the cut y r - ’ > y 
in the J case than in the D and C ones. 

We do not push any further similar comparisons involving the J scheme. We have 
retained it in Fig. fj] because it has been used extensively at LEP1 in the past years. However, 
such algorithm has been proven to be more sensitive to hadronisation effects than the D 


and C algorithms, especially at low y values, where the six-jet fraction onsets |8f 
Furthermore, the Jade rates do not show the usual Sudakov exponentiation of multiple soft- 
gluon emission |38|], despite having an expansion of the form a” ln m (l/y) at small values 
of the resolution parameter, contrary to the case of the Durham and Cambridge schemes 







(see [^3] and references therein), and generally have a scale dependence of the fixed-order 
perturbative rates (e.g., the three- |85|] and four-jet [^(J fractions) larger than that of the 
D and C schemes. Thus, given the mentioned flaws, in the following we will give total and 
differential rates only for the last two case^]. 

As for these, one should notice that the larger rate for the C scheme as compared to the 
D one is a direct consequence of the ‘soft freezing’ procedure of resolved jets described in 


Ref. |L|. The step of eliminating from the sequence of clustering the less energetic one in 
a resolved pair of particles (i.e., with > y ), implemented in the C algorithm, tends to 
enhance the final jet multiplicity of the original D scheme. In fact, the procedure prevents 
the attraction of the remaining particles (at large angle) into new unresolved pairs, whose 
momenta would then be merged together, producing a lower number of final jets. For 
example, at y = (0.001)[0.005]{0.010}, relative differences between the two algorithms are 
of the order of (10) [14] {16}%. 

The jet fraction / 6 (j/) and cross section <r 6 (y) corresponding to the rates in Fig. £| are given 
in Tab. [5]. for a representative selection of y’s. The value of a s adopted is 0.120, whereas 
that used for at is 39.86 nb. Given the total luminosity collected at LEP1 in the 1989-1995 
years (more than 10' hadronic events have been recorded by the four Collaborations), the 
six-jet event rate is comfortably measurable. For example, for a luminosity of, say, 100 
pb -1 per experiment, perturbative QCD predicts at LO some 58,000 original six-parton 
events recognised as six-jet ones, for y = 0.001 in the C scheme (a number that decreases 
to approximately 52,000 if the D algorithm is adopted instead). 

The rates given in Figs. f| and Tab. are LO results. Next-to-leading order (NLO) 
corrections proportional to O(a^) are expected to be large. In fact, the size of higher-order 
(HO) contributions generally increases with the power of a s and with the number of particles 
in the final state or, in other terms, with their possible permutations (i.e., with the number 
of different possible attachments of both the additional real and virtual particles to those 
appearing in lowest order). As already mentioned, the highest-order corrections calculated 
to date in e + e _ —> n-jet annihilations are the NLO ones to the four-jet rate, that is, terms 
proportional to ctf 


14| . The total four-jet cross section at NLO was found to be larger 
than that obtained at LO in by a K-factor of 1.5 or more, depending on the value of 
y implemented during the clustering procedure and the algorithm used as well. Therefore, 
we should expect that the rates given in Tab. § underestimate the six-jet rates by at least 
a similar factor. Indeed, assuming a K cs 2 correction, our results are in a satisfactory 
agreement with both data and phenomenological Monte Carlo (MG) programs [|87| . 

10 Apart from one example, when we will need to compute a cross section using the J scheme for purpose 
of comparison against results published in literature. 









y 

fe(y) 

^e(y) (pb) 

0.001 

(1.31) [1.47] x 10 - 2 

(523.91) [584.19] 

0.002 

(3.16) [3.65] x 10 “ 3 

(125.95) [145.44] 

0.003 

(1.17) [1.37] x 10 ~ 3 

(46.81) [54.72] 

0.004 

(5.29) [6.01] x 10 “ 4 

(21.09) [23.96] 

0.005 

(2.63) [3.05] x 10 ” 4 

(10.49) [12.17] 

0.006 

(1.38) [1.59] x 10 " 4 

(5.49) [6.35] 

0.007 

(7.90) [8.85] x 10 " 5 

(3.15) [3.53] 

0.008 

(4.80) [5.20] x 10 " 5 

(1.91) [2.07] 

0.009 

(2.77)[3.08] x 10 ~ 5 

(1.10)[1.23] 

0.010 

(1.64) [1.93] x 10 " 5 

(0.65) [0.77] 

LEP1 


Table 5: Jet fraction and cross section rates for e + e _ —> 6 jets at LEP1, for several values 
of y in the Durham and Cambridge schemes, in round and squared brackets, respectively. 
The numerical errors do not affect the significant digits shown. 


In this respect the MCs J88|-||92|] largely exploited in experimental analyses apparently 
perform better than the fixed-order perturbative calculations, especially at low ^/-values, 
where the latter need to be supported by the additional contribution of perturbation series 
involving powers of log(l/y) resummed to all orders (for the case of the D and C scheme, 
see Refs. [[I9|, |40], ^3|) in order to fit the same data. However, as such MCs generally 
implement only the infrared (i.e., soft and collinear) dynamics of quarks and gluons in the 
standard ‘parton-shower (PS) + 0(a s ) Matrix Element (ME)’ modelling, in many cases their 
description of the large y-behaviour and/or that of the interactions of secondary ‘branching 
products’ is (or should be expected to be) no longer adequate. In fact, this has been shown 
to be the case, e.g., for some typical angular quantities of four-jet events E3- H 


As 


discussed in Ref. , these differences are presumably due to a poor description of the spin 
correlations among the various partons implemented in the MCs: in particular, between the 
primary (i.e., from the 7 *, Z^-splitting) and secondary (i.e., from the gluon splitting into 
gluons or quarks) pairs. In contrast, once 0(a 2 ) MEs are inserted and properly matched 
to the hadronisation stage, e.g., using the JETSET string fragmentation model [|0(J (see 
also Ref. |)5]]), then the agreement is recovered! 11 ]. If one considers that four-jet events 


11 In this context, one should however mention that ME models with ‘added-on’ fragmentation, i.e., l O(a^) 
+ hadronisation’ cannot be reliably extrapolated from one energy to another, as the fragmentation tuning 
is energy dependent. To obviate this, a special l O(al) + PS + cluster hadronisation’ version of HERWIG is 






























represent only the next-to-lowest-order QCD interactions in e + e _ scatterings, then it is 
not unreasonable to argue that further complications might well arise as the final state 
considered gets more and more sophisticated, such as in six-jet events. 

In order to emphasise this point, we present in Figs. the differential distributions of 
the six-jet rate (along with the five- and four-jet ones) in some typical shape variables used 
in QCD studies in electron-positron annihilations. They are the thrust T 0 , the oblateness 
, the C-parameter C [Q, the heavy Mh , light Ml and difference Mp jet masses ||100 


O 


the sphericity S, aplanarity A and planarity P [|101|1 , the minor m and major M values |j 8 


and the D-parameter [[hj . (Their allure in testing the underlying QCD dynamics has been 
well illustrated in, e.g., Ref. [0].) In Figs. |5|-§|the algorithm used to produce the plots was 
the Cambridge one, with y = 0.001. 

As compared to the same distributions in the case of lower jet multiplicity events (such 
as, e.g., four-jet production), those in Figs. 0-0 for the six- (more markedly) and five-jet (in 
some instances) rates are much more ‘central’. That is, the peaks of the distributions are 
generally further away from the infrared regions, those which dominate the phenomenology 
of rz.-jet events, for n = 3 and 4. (The only exception is the C-parameter, which is indeed 
expected to play a role apart since there exists a singularity at the boundary of the three- 
parton phase space C = 2/3, see Ref. [ |102|| , affecting the four-jet distribution but not the 
five- and six-jet ones.) This is a clear indication that other dynamics (other than the soft 
and collincar emission of gluons off quark lines) has a perceptible influence on differential 
jet rates, as the multiplicity of the latter increases. In other terms, one should expect 
the above mentioned spin correlations among the partons to modify quite significantly the 
infrared behaviour dictated by the propagators when the energy of the particles decreases, 
as happens in the six-parton case and particularly for the secondary and tertiary branching 
products of the g —> qq, g —> gg and g —» ggg splittings (see Fig. [[]). 

To further stress this, we also plot in Figs. [o| the same shape variables distributions 
as obtained at hadron level for the six-jet (i.e., n — 6 ) rate by using the HERWIG event 
generator, with identical choice of jet algorithm and resolution as above. As a matter of 
fact, the MC spectra depart in most cases quite substantially from those obtained by using 
the exact O(a^) MEs, notably, they ‘point’ towards the singularities. Such behaviour is 
particular evident for thrust and sphericity, two variables to which we will come back below. 

Under these circumstances, we believe the availability of exact perturbative calculations 
of six-jet events to be essential in order to model correctly the HO parton dynamics for 
both future tests of QCD (like those outlined here) and QCD background studies (as we 
will show later on for the case of tt —> 6 -jet production and decay at NLC). 


now in preparation [p6|. 


























The total hadronic cross section falls drastically when increasing the CM energy from 
the LEP1 to the LEP2 values, by more than three orders of magnitude, and so does the 
six-jet rate. Nonetheless, six-jet fractions have been measured also during the 1995-1996 
runs at the CM energies of 130-136, 161 and 172 GeV, when a total luminosity of around 27 
pb -1 was collected, and are currently under further investigation (at 183 GeV and higher). 


Results can be found for the lower energies in Ref. ||103|| . Within a ‘typical’ K-factor of 2 
(which tentatively quantifies the ratio between the NLO and the LO six-jet rates also at 
LEP2) the values of as reconstructed from our rates are always well compatible with 
those produced by the MCs used by e.g., the ALEPH Collaboration ||104|] , (Notice that this 
is also true at y/s = 161 GeV, where a six-jet ‘excess’ has apparently been observed ||104j| ). 
For reference, we present the jet rates at yfs = 180 GeV in Fig. [7] in the form of the total 
cross section. The relative differences between the two D and C algorithms are similar to 
those already encountered at LEP1. 

Quite apart from its relevance in QCD jet fractions, the six-jet rate at LEP2 can in 
principle represent a source of background to the so called ‘colour reconnected’ W + W~ 
hadronic decays ||105|| . The problem goes as follows. Colour exchange (both at perturbative 
and non-perturbative level) can take place between the decay products of the two W^’s, 
thus destroying the notion of two separate W ± decays and hence casting doubts on our 
abilities to reconstruct accurately the mass of the resonances. In perturbative terms, such 
effect is due to the presence of interference terms in the total cross section for W + W~ —> 
multi-parton production and decay, which account for the exchange of real or virtual partons 
emitted in the decay of one of the IE ±, s and absorbed by the partons generated by the other. 
Simple rules of colour conservation show that such colour effects cannot occur at first order 
in a SJ but only in higher orders (i.e., at least two gluons need to be exchanged). At the 
non-perturbative level, a similar phenomenon takes place when the colour strings generated 
among the partons from the Erst W ± decay overlap with those from the second one, break 
down and form new strings which may join quarks from different shower cascades (e.g., 
in the language of the Lund fragmentation model |U06|| ). Of the two contributions, the 
non-perturbative one is expected to be dominant ||107|| . The consequence is similar at both 
levels: namely, particles are generated that cannot be unambiguously assigned to one W ± 
decay or another. 

It is of decisive importance to quantify accurately the above effects, especially if one 
wants to achieve a high precision measurement of the mass M w of the W ± boson (a target 
accuracy of ±(40 — 50) MeV or less is expected, see, e.g., Ref. [ |108|| ). In fact, at present, 
their dynamics is relatively unknown. On the one hand, a complete O(a^) perturbative 
calculation of W + W~ production and decay does not exist to date, as only the tree-level 
part of e + e _ —> W + W~ —> 6 jets has been computed so far [[53]. On the other hand, in the 






























non-perturbative modelling, many more additional parameters are needed to describe the 
reconnection effects than the well tested mechanisms of the quark and gluon fragmentation, 
so that there are many phenomenological programs now on the market, indeed predicting 
rather different sizes for the above colour effects m- 


Whereas the treatment of non-perturbative colour rearrangement effects is beyond the 
scope of the paper, we come back to the perturbative ones. In particular, we would like 


to recall that Ref. |5£| gave an estimate of the size of the latter at LO, in the form of the 
part of the e + e _ — > W + W~ —> 6 -jet cross section due to the mentioned interference terms. 
The effect was found to be quite small: e.g., at the level of 4 x 1CT 3 pb at y — 0.001 
for the Jade algorithm, with y/s = 175 GeV and using a s = 0.115. The cumulative rates 
for processes ( 0 ) _ ( 0 ) f° r the same jet clustering scheme and the identical combination of 
y, y/s and a s amount to approximately 2 pb, well above the colour reconnected rates in 
W + W~ —» 6 -jet events. Since in the latter case colour exchange proceeds via gluons, in 
Ref. |Q the corresponding energy spectra were studied, in order to give some hints about 
the sort of energy at which perturbative colour reconnection takes place. The softest gluon 
energy distribution was found to peak at about 8 GeV, whereas that of the most energetic 
one lies at approximately 20 GeV. To a good degree of approximation, the two gluon energy 
spectra coincide with those of the two least energetic particles in W + W~ —» 6 -jet decays. 
In Fig. | 8 |, we have studied the energy spectra of the partons produced via reactions (|1])-(|3|). 
Coincidentally enough, the fifth and the fourth most energetic ones show their maxima 
around 8 and 20 GeV, respectively^. That is, they behave on the same footing as the 
gluon jets in e + e _ —> W + W~ —> 6 -jet events, which are responsible for perturbative colour 
rearrangement. Though both the studies in Ref. |58|] and those presented here need to 
be eventually supported by loop-calculations as well, it is clear that six-jet events from 
QCD represent a considerable source of irreducible background to hadronic W + W~ events 
perturbatively reconnected. 

As a matter of fact, one might consider imposing a cut such as Mjj ~ M Wl that is, asking 
that one di-jet invariant mass (see below, eq. ([29]) ) is close to the value of the W ± mass, as 
one might naively expect to occur naturally in W + W~ —> 6 -jet events. However, this could 
well be counter-productive, for two reasons. Firstly, the colour reconnected contribution to 
e + e~ —> W + W~ —> 6 -jet events is defined as the sum of two terms, namely, the interference 
between the set of diagrams labelled (Al) and (A2) plus that between (A3) and (A4) 
in Fig. EQ (Note that the cross term between (Al)©(A2) and (A3)©(A4) is identically 

12 Further notice that also in events of the type ® (§ the fifth and fourth most energetic partons are 
in most cases gluons, because of the dominance in the total QCD cross section of two-quarks-four-gluon 
events, with the two fermion being the most energetic ones. 

13 We would like to thank the authors of Ref. 158 for their kind permission of exploiting here one of the 













zero). Neither of the two terms has a distinctive peak at Mjj « Mw, in line with the notion 
of colour exchange taking place between the two W ± decays. Secondly, the invariant mass of 
the di-jet combination formed by pairing the first and third most energetic jets from six-jet 
QCD events peaks not far from the W ± mass, i.e., at ~ 76 GeV (less than two IV from 
M\y) and is rather broad in the region Mjj ~ Mw, as can be appreciated in Fig. [10] (here, 
our default algorithm, resolution parameter and CM energy for LEP2 have been used again, 
as the shape of the distribution is not significantly affected by either of those, provided 
y/s 2 Mw for the energy). Further notice that other QCD di-jet mass distributions 

populate the region around the Mw resonance, such as the combinations ij = 12,14 and 23 
(the other spectra are much softer, as they concentrate at low invariant masses). Evidently 
then, under these circumstances, our understanding of perturbative colour interconnection 
effects in W + W~ —> 6-jet decays is bound to our ability in suppressing the corresponding 
QCD background. This operation can however be attempted only in the context of a detailed 
Monte Carlo simulation, including parton shower, hadronisation and detector effects, which 
is beyond our capabilities. 

The six-jet event rate that will be collected at NLC can be rather large, despite the cross 
section being more than a factor 10 4 smaller than at LEP1 (e.g., at y/s = 500 GeV). This 
is due to the large yearly luminosity expected at this machine, around 10 fb _1 . For such a 
value and assuming a standard evolution of the coupling constant with the increasing energy 
(that is, no non-Standard Model thresholds occur up to 500 GeV), at the minimum of the 
y -values considered here (i.e., y = 0.001), one should expect some 220 events per annum by 
adopting the D scheme and about 13% more if one adopts the C one. However, these rates 
decrease rapidly as the resolution parameter gets larger, by a factor of 50 or so at y = 0.005 
and of approximately 800 at y = 0.01. 

It is also interesting to look at the composition of the total rates in terms of the three 
processes ©-(!)• Whereas this is probably of little concern at LEP1 and LEP2, the capabil¬ 
ity of the detectors of distinguishing between jets due to quarks (and among these, bottom 
flavours in particular: e.g., in tagging top and/or Higgs decays) and gluons, is essential 
at NLC, in order to perform dedicated searches for old and new particles. The separate 
behaviours of the reactions in ©-(I!) can be appreciated in Fig. [□], e.g., for the case of the 
C scheme, in terms of total cross sections. The rates for the D algorithm follow the same 
pattern. 

The tagging of 6-quarks produced, e.g., in the hadronic top decay chain tt —> 
bbW + W~ —> bbjjjj, is of decisive importance, non only to suppress the QCD background 
(mainly proceeding via gluon events, see Fig. B but also to reduce the combinatorics in¬ 


figures of their paper. 



volved in selecting three jets out of six, in case no flavour identification is exploited (twenty 
combinations in total). If 6 -quark jets can be recognised (via their semileptonic decays, via 
/i-vertex devices, etc.), then only six three-jet combinations bjj (or bjj ) need to be selected, 
one of which is such that ~ m t . If is the efficiency of tagging one single 6 -jet 

within a certain fiducial volume, and two taggings are required, than e 2 is the overall rate 
for identifying events containing exactly two 6 -quarks (we assume no correlations between 
the two 6 -taggings and that the charge of the heavy quark is recognised). This is the case 
which applies to the subprocesses e + e” —» bbgggg and e + e” —» bbqqq'q' (with q, q' % &)Q 
The probability of misidentifying six-jet events produced via four- 6 -quark reactions (i.e., 
e + e _ —> bbbbgg or e + e _ —> bbbbqq (with q % 6 )) as bbjjjj events is then [2ef,(l — e;,)] 2 . If 
there are six 6 -quarks in the event, as in e + e _ —> 666666 , then the corresponding number 
is [3eb(l — Cfe) 2 ] 3 - The total rates as a function of y in the C scheme for six-parton process 
involving 6 -quarks are given in Fig. B for the case in which only u and d massless contri¬ 
butions are considered (for c and s, respectively, one gets very similar results, even when 
the masses, e.g., m c = 1.35 GeV and m s = 0.30 GeV, are retained in the calculation). 

If one multiplies the rates given in Fig. [L2]by the overall 6 -tagging efficiencies mentioned 
in the previous paragraph, then one realises that six-jet events produced via four- and six- 
6 -quarks subprocesses yield negligible rates at NLC (assuming, e.g., e& = 70%, the current 
LEP value). Events involving only two 6 -quarks can instead produce a detectable number of 
bbjjjj events, especially at low y. For example, for y = 0.001 in the C scheme one gets some 
15 QCD background events in bbjjjj samples (including the 6 -tagging efficiency el), after a 
luminosity of ten inverse femtobarns has been collected. Though this is not a large number 
per se, considering that e + e _ —> tt —» 6 jets yields about 800 doubly 6 -tagged events at the 
same t/-value in the above scheme, one should recall that very high precision measurements 
of top parameters (such as m t , T t , the branching ratios, etc.) are expected to be carried out 
at NLC m- Under these circumstances, even small backgrounds need to be quantified 


accurately. 

Whereas mass effects (markedly, of 6-quarks) are always negligible (at LEP1, LEP2 and 
also NLC) if the sum over all flavour combinations is performed (that is, when no heavy 
flavour tagging is exploited), they can in principle be important in 6-jet samples, as in 
those selected for top studies at NLC. As a matter of fact, in the individual contributions 
to the six-jet rate mass corrections can be quite sizable at a/s = 500 GeV, as shown in 

14 For simplicity, we do not consider here the case of c-quark misidentified as 6-quarks. We can presumably 
assume that by the time NLC will be operative the efficiency in distinguishing between c- and 6-quark jets 
will be much higher than at present. The rejection factor against light quark jets (produced by u-, d- and 
s-quarks) as well as against gluon jets is very high already at present, so that we do not need to worry 
about their contamination in the present analysis. 




Fig. Q. at the level of ten percent or more. These corrections are more visible when more 
6 -quarks are produced in the final state. However, as previously mentioned, these multi- 
6 -quark subprocesses are unlikely to be important in top analyses. In practise, the largest 
part of the bbjjjj sample produced from QCD comes from two-quark-four-gluon events, 
where 6 -mass effects are naturally smaller because only two massive particles are involved, 
over a six-particle phase space and at rather high energy (y/s ^ 2 mj). In the end, 6 -mass 
effects on the complete bbjjjj produced by QCD amount to only 4% at the NLC, for typical 
values of the resolution parameter in the six-jet region, y « 0.001. We have also verified 
that the differential distributions (e.g., energy and angles) suffer imperceptibly from mass 
corrections, so is the case of the shape variable spectra. Therefore, in first instance, it is 
safe to assume the value mj = 0 in numerical simulations of QCD six-jet events. We will 
do so in the remainder of the paper. 

The concern about possible background effects at the NLC due to six-jet events from 
pure QCD comes about if one further considers that they may naturally survive some of the 
top signal selection criteria. These can in fact rely on the analysis of shape variable spectra 
of six-jet samples, in particular, the thrust and sphericity distributions, and since the latter 
(as already mentioned while commenting Figs. §0) might well be poorly described by the 
parton shower approach exploited in the phenomenological MC programs, it is of decisive 
importance that also the non-infrared dynamics of events of the type (|])-(§) is carefully 
studied. 

It is well known that the large value of the top mass leads to rather spherical signal events. 
Therefore, shape variables such us the two mentioned above represent useful means to 
disentangle e + e~ — > tt events. For example, a selection strategy that does not exploit neither 
lepton identification nor the tagging of 6 -jets was outlined in Ref. ||111|| . This approach has 
a twofold attractiveness: firstly, its simplicity; secondly, it is a rather general procedure, 
that can be equally applied both at and above the threshold y/s ~ 2 m t . The requirements 
are a large particle multiplicity, a high number of jets (at least five, eventually forced to six 
by a clustering scheme) and a rather low thrust (typically, below 0.85). Jets are selected 
according to a jet clustering algorithm (in our case, the Cambridge one with y = 0.001, for 
sake of illustration). 

Clearly, the six-jet events meet the two first criteria. As for the thrust and sphericity 
distributions, these are shown in Fig. [L3]. From there it is evident the overlapping of the two 
distributions {tt and QCD events) in the region where the signal events are mostly located, 
for both thrust and sphericity. A cut such as that advocated in Ref. ||111|| , i.e., T < 0.85, 
would then be rather ineffective against six-jet events from QCD. Once again, we believe 
that a correct modelling of the phenomenology of the latter based on a ME description will 
be essential to quantify correctly the effects of the QCD background, especially in view of 









the mentioned high precision measurements of top parameters. 

Tests of multiple electroweak self-couplings of gauge-bosons (as well as searches for 
new resonances) at NLC will often need to rely on the mass reconstruction of multi-jet 
systems (particularly di-jet ones). Therefore, it is instructive to look at the invariant mass 
distributions which will be produced by all the possible two-parton combinations ij in six- 
jet events from QCD at 0(a ^) (with i — 1, ...5 and j — i + 1, ...6). As usual in multi-jet 
analyses, we first order the jets in energy, so that E\ > E 2 > ... > E 5 > E§. Then, we 
construct the quantities (the equality strictly holds only for massless partons ij, as is the 
case here) 


Ml 


nriij = 


^i 3 


2EiEj(l — cos Oij ) 


(29) 


where represents the Lorentz invariant (squared) mass. These fifteen quantities are 
shown in Fig. [14| for the C scheme at y — 0.001, their shape being similar for the D 
algorithm. We found it convenient to plot the ‘reduced’ invariant masses rriij rather than 
the actual ones M i3 , as energies and angles ‘scale’ with the CM energy in such a way that the 
shape of the distributions is largely unaffected by changes of the value of a/s in the energy 


range relevant to NLC. Therefore, from Fig. |14| one should then be able to reconstruct rather 
accurately the Lorentz invariant mass distributions for a given CM energy a /s by exploiting 
eq. (|25|). Note that the reduced mass spectra are similar in all three reactions (jf])-([3p, their 
integral being however rescaled according to the numbers given in Fig. 0 To allow for an 
easy conversion of the differential cross sections into numbers of events in a certain mass 
range, the spectra in Fig. [F| sum to the total cross section of processes (|1])-(|3|) NLC. 

It is interesting to notice in Fig. [L4|the ‘resonant’ behaviour of some of the distributions. 
This is particularly true for those involving the most energetic of all the partons. Using 
the definition (^9|) they translate into peak-like structures at the true invariant mass values 
Mij ~ (250)[212]{177} GeV, for the combinations ij = (12)[13]{14}, and, possibly, M 15 ~ 
150 GeV as well. In all the other cases the spectra are generally softer and do not show 
any distinctive kinematic feature. Once again, it is important that all such behaviours 
are correctly implemented in the simulation programs that will be adopted by the NLC 
experiments, so that it will be possible to recognise and eventually subtract the six-jet 
QCD background. 

A further general comment that is in order for six-parton QCD events produced at NLC 
is that the much lower value of a s at those high energies (as compared to that at LEP1 
and/or LEP2) in principle implies a reduced importance of the uncalculated HO strong 
corrections, so that the theoretical uncertainty on the rates of processes ©-( 0 ) will be 
more under control at the future collider than it was/is at the energy scales of the two LEP 
phases. 





So far we not included effects of ISR in any of our calculations for six-jet production. 


However, it is well known j53] that photon bremsstrahlung generated by the incoming 
electron and positron beams can be quite sizable at NLC. In contrast, at LEP1 and LEP2 
such radiation is suppressed by the width of the Z and W + W~ resonances, respectively^. 
In order to implement ISR we resort to the mentioned ESF approach. In particular, we use 
the expressions given in Ref. [|1 14|| , through the order a 2 em . We plot the production rates for 
the sum of the processes fl!])-© in presence of ISR in Fig. [15], as a function of the resolution 
parameter y in the C scheme. (Effects turn out to be similar in the D scheme.) They are 
compared to those obtained without ISR (that is, the sum of the rates in Fig. [TT]) . No cut is 
here applied to the soft and collincar photons from the initial state. We see that the curve 
corresponding to the processes (0)-© convoluted with ISR lies above the lowest-order one. 
This is rather intuitive, as it is well known that the radiation of photons from the incoming 
electron and positrons tends to lower the ‘effective’ CM energy of the collision ||1 13 


Since 


processes (0)-© proceed via s-channel diagrams, the expected overall effect would be an 
enhancement of the total rates. In fact, this is what can be seen in Fig. |T^. At y = 0.001, 
the difference is around 25% and this tends to increase as y gets larger. However, we have 
verified that the shape of the differential distributions (such as those in shape variables, see 
Fig. |T3|. and in invariant mass, see Fig. H) suffer little from ISR effects. 

As announced in the Introduction, we also devote some space to analyse photoproduction 
in 5jet + ly, 4jet + 2y, 3jet + 3y and 2jet + 4y events. For this study, we reconsider the 
case of LEP1 energies, where the impact of ISR is much less relevant than at LEP2 and 
NLC (recall that we did not implement the resummed logarithmic part of it in photonic 
events). I 11 the first two cases we consider both the dominant contributions involving only 
two quarks and the suppressed ones involving four, whereas in the last two cases only two 
quarks can participate in the event. 

Final states including photons are actively studied experimentally in order to determine 
the electroweak couplings of the quarks and as a mean to search for new phenomena ||115| . 


Furthermore, such particles are the only ones which can be directly revealed and couple to 
the quarks at each stage of the perturbative parton shower (see Fig. []]), so that the detection 
of hard isolated photons (or ‘direct photons’) is extremely useful for a detailed understanding 
of the underlying partonic picture as well as of the QCD evolution from large to small energy 
scales. (Conversely, soft photons lying close to hadronic tracks can give information about 
the parton fragmentation mechanisms, though the corresponding rates are not calculable 


15 In addition, peculiar features at NLC is the presence of Linac energy spread and beamsstrahlung |112| . 
However, their effects are negligible, as compared to those produced by ISR, at least for the ‘narrow’ D-D 
and TESLA collider designs [ 113 ]. For this reason, we do not consider them in our analysis. They would 


however be straightforward to insert. 



















perturbatively.) In particular, the more jets are separated in association with photons, the 
more likely will be that photons were radiated by quarks in later stages of such evolution. 
In e + e _ —> m jet + n photon events (with m + n = 6 and m > 2), 7 -radiation always takes 
place from the primary quark-antiquark pair (i.e., that produced in the 7 , Z-splitting) in 
the case of the dominant two-quark events. Thus, for the cases m = 4(5) and n = 2(1), 
such contributions should be regarded as a source of background in the mentioned analyses, 
with the four-quark contributions being the interesting ones, as in these cases photons can 
also be radiated by secondary branching products. 

In this study, we limit ourselves to giving some total rates. The cross sections for isolated 
photon events as a function of the y resolution parameter, applied to both quarks/gluons and 
photons on the same footing (in the spirit of the so-called ‘democratic approach’ | T16|| ), in the 
C scheme, are presented in Fig. [16], at LEP1. The mentioned cuts in transverse momentum 
and polar angle of the photons have been enforced, in order to screen the singularities due to 
infrared (i.e., soft and collinear) EM emission from electrons and positrons. The hierarchy in 
the cross sections of Fig. [TO] is dictated mainly by the perturbative order 0(a^ 2 a™) for the 
hard production of m jets and n photons. At LEP1, given, e.g., some 400 inverse picobarn 
of total accumulated luminosity, one should expect that only the 5jet + I 7 and 4jet + 2y 
events can yield an observable rate, with the four-quark contributions being roughly one 
order of magnitude smaller than the two-quark ones in both cases. 

4 Summary 

In summary, we have studied at the partonic level the tree processes e + e~ —> qqgggg , 
e + e _ —* qqq'q'gg and e + e - —> qqq'q'q"q" (for massless and massive quarks) at leading-order 
in perturbative QCD, by computing their exact matrix elements. The latter task has been 
accomplished thanks to the use of spinor techniques and colour decomposition methods. 
Both the formulae for the Feynman amplitudes and the relevant QCD (i.e., for Nq = 3) 
colour factors have been given explicitly. Several phenomenological applications at present 
and future colliders have been discussed. Two FORTRAN codes have been produced and cross 
checked with each other on several different platforms: VMS, OSF, HP. The fastest of the 
two requires about 3 x 10 ~ 3 CPU seconds to evaluate a single event on an alpha-station 
DEC 3000 - M300X, which makes it suitable for high statistics Monte Carlo simulations. 
Further optimisations are, however, in progress. 



5 Acknowledgements 


We acknowledge the UK PPARC support. We thank the Theoretical Physics Groups in 
Fermilab and Lund for their kind hospitality while part of this work was carried out. This 
research work was supported in part by the Italian Institute of Culture ‘C.M. Lerici’ under 
the grant Prot. I/Bl 690, 1997. Finally, we are grateful to David J. Miller for carefully 
reading the manuscript and for several useful discussions and comments. 


References 

[1] R.K. Ellis, W.J. Stirling and B.R. Webber, “QCD and Collider Physics” (Cambridge 
University Press, 1996). 

[2] M.L. Mangano and S.J. Parke, Phys. Rep. 200 (1991) 301. 

[3] R. Gastmans and T. T. Wu, “The ubiquitous photon” (Addison Wesley, 1990). 

[4] F. Bloch and A. Nordsieck, Phys. Rev. 52 (1937) 54. 

[5] T. Kinoshita, J. Math. Phys. 3 (1962) 650; 

T. D. Lee and M. Nauenberg, Phys. Rev. 133 (1964) 1549. 

[6] G. Sterman, Phys. Rev. D17 (1978) 2789. 

[7] See, e.g.: 

L. M. Brown and R. P. Feynman, Phys. Rev. 85 (1952) 231; 

D. B. Melrose, R Nuovo Cimento 40A (1965) 181; 

A.C.T. Wu, Mat. Fys. Medd. Dan. Vid. Selsk. 33 No. 3 (1961) 72; 

G. Passarino and M. Veltman, Nucl. Phys. B160 (1979) 151; 

G. ’t Hooft and M. Veltman, Nucl. Phys. B153 (1979) 365; 

W. van Neerven and J.A.M. Vermaseren, Phys. Lett. B137 (1984) 241; 

G.J. Van Oldenborgh and J.A.M. Vermaseren, Z. Phys. C46 (1990) 425; 

A. Denner, LI. Nierste and R. Scharf, Nucl. Phys. B367 (1991) 637; 

A. Denner, Fort. Phys. 41 (1993) 307; 

U. Nierste, D. Muller and M. Bohm, Z. Phys. C57 (1993) 605. 

[8] See, e.g.: 

Z. Bern, L. Dixon and D.A. Kosower, Nucl. Phys. B425 (1994) 217. Nucl. Phys. Proc. 
Suppl. 51C (1996) 243; Ann. Rev. Nucl. Part. Sci. 46 (1996) 109; Nucl. Phys. B, Proc. 
Suppl. 39BC (1995) 146; 



Z. Bern, L. Dixon, D.C. Dunbar and D.A. Kosower, Nucl. Phys. B435 (1995) 59; talk 
presented at ‘Continuous Advances in QCD’, Minneapolis, Feb. 18-20, 1994, preprint 
SLAC-PUB-6490, May 1994, hep-ph/9405248; talk presented at ‘DIS ’97’, Chicago, 
IL, April 14-18, 1997, preprint Saclay-SPhT/T97-62, June 1997, hep-ph/9706447. 

[9] N.K. Falck, D. Graudenz and G. Kramer, Phys. Lett. B220 (1989) 299; Nucl. Phys. 
B328 (1989) 317; 

K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B313 (1989) 560; 

F.A. Berends, W.T. Giele and H. Kuijf, Nucl. Phys. B321 (1989) 39. 

[10] R.K. Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B178 (1981) 421. 

[11] See, e.g.: 

F. Mandl and G. Shaw, “Quantum Field Theory” (John Wiley & Sons, 1993). 

[12] S. Moretti, Phys. Lett. B420 (1998) 367. 

[13] E.W.N. Glover and D.J. Miller, Phys. Lett. B396 (1997) 257; 

J.M. Campbell, E.W.N. Glover and D.J. Miller, Phys. Lett. B409 (1997) 503. 

[14] Z. Bern, L. Dixon, D.A. Kosower and S. Weinzier, Nucl. Phys. B489 (1997) 3; 

Z. Bern, L. Dixon, D.A. Kosower, Nucl. Phys. B513 (1998) 3; 

A. Signer, talk presented at XXXIInd Rencontres de Moriond: QCD and High-Energy 
Hadronic Riteractions, Les Arcs, France, 22-29 March 1997, preprint SLAC-PUB-7490, 
May 1997, hep-ph/9705218; 

L. Dixon and A. Signer, Phys. Rev. D56 (1997) 4031. 

[15] W. Bernreuther, A. Brandenburg and P. Uwer, Phys. Rev. Lett. 79 (1997) 189; 

A. Brandenburg and P. Uwer, Nucl. Phys. B515 (1998) 279; 

G. Rodrigo, Nucl. Phys. Proc. Suppl. 54A (1997) 60; 

G. Rodrigo, A. Santamaria and M. Bilenky, Phys. Rev. Lett. 79 (1997) 193; 

P. Nason and C. Oleari, Nucl. Phys. B521 (1998) 237. 

[16] B.L. Ioffe, Phys. Lett. B78 (1978) 277. 

[17] G. Branco, K.H. Streng and H.P. Nilles, Phys. Lett. B85 (1979) 269. 

[18] A. Ali, J.G. Korner, Z. Kunszt, E. Pietarinen, G. Kramer, G. Schierholz and J. Willrodt, 
Nucl. Phys. B167 (1980) 454; 

J.G. Korner, G. Schierholz and J. Willrodt, Nucl. Phys. B185 (1981) 365. 


[19] K.J.F. Gaemers and J.A.M. Vernraseren, Z. Phys. C7 (1980) 81. 



[20] O. Nachtmann and A. Reiter, Z. Phys. C16 (1982) 45. 

[21] D. Danckaert, P. de Causmaeker, R. Gastmans, W. Troosts and T.-T. Wu, Phys. Lett. 
B114 (1982) 203. 

[22] L. Clavclli and G.V. Gehlcn, Phys. Rev. D27 (1983) 1495. 

[23] F.A. Berends, P.H. Daverveldt and R. Kleiss, Phys. Lett. B148 (1984) 489; Nucl. Phys. 
B253 (1985) 441. 

[24] E. Laermann, K.H. Streng and P.M. Zerwas, Z. Phys. C3 (1980) 289. 

[25] J. Jersak, E. Laermann and P.M. Zerwas, Phys. Rev. D25 (1982) 1218. 

[26] E.W.N. Glover, R. Kleiss and J.J. van der Bij, Z. Phys. C47 (1990) 435. 

[27] E.N. Argyres, C.G. Papadopoulos and S.D.P. Vlassopulos, Phys. Lett. B237 (1990) 
581. 


[28] G. Kramer and B. Lampe, Fort. Phys. 37 (1989) 161. 

[29] A. Ballestrero, E. Maina and S. Moretti, Phys. Lett. B294 (1992) 425; Nucl. Phys. 
B415 (1994) 265; Proceedings of the XXIXth Rencontres de Moriond: QCD and High 
Energy Hadronic Interactions, Meribel, Savoie, France, March 1994, ed. by J. Tran 
Thanh Van, ed. Frontieres, Gif-snr-Yvette, 1994, 367. 


[30] A. Ballestrero and E. Maina, Phys. Lett. B323 (1994) 53. 


[31] W.T. Giclc and E.W.N. Glover, Phys. Rev. D46 (1992) 1980; in Ref. fy], Vol. 2, page 
163. 


[32] Z. Kunszt and P. Nason, Proceedings of the Workshop L Z Physics at LEP P, ed. by 
G. Altarclli, R. Kleiss and C. Verzegnassi, CERN Yellow Report No. 89-08, Vol. 1, page 
373; in Ref. [ET1, Vol. 2, page 163. 


[33] S. Catani and M.H. Seymour, in Ref. [511] , Vol. 2, page 163. 

[34] Z. Nagy and Z. Trocsanyi, Phys. Rev. Lett. 79 (1997) 3604; Phys. Lett. B414 (1997) 
187; talk presented at the Euroconference on Quantum Chromodynamics: QCD 97: 
25th Anniversary of QCD, Montpellier, France, 3-9 July 1997, preprint C97-07-03.1, 
August 1997, hep-ph/9708344. 


[35] A. Signer, Comp. Phys. Commun. 106 (1997) 125. 











[36] J.M. Campbell, M. Cullen and E.W.N. Glover, in preparation. 

[37] S. Catani, G. Turnock, B.R. Webber and L. Trentadue, preprint Cavendish-HEP-90/16, 
August 1990 (unpublished); Nucl. Phys. B407 (1993) 3; 

G. Dissertori and M. Schmelling, Phys. Lett. B361 (1995) 167. 

[38] N. Brown and W.J. Stirling, Phys. Lett. B252 (1990) 657. 

[39] S. Catani, Yu.L. Dokshitzer, M. Olsson, G. Turnock and B.R. Webber, Phys. Lett. 
B269 (1991) 432. 

[40] N. Brown and W.J. Stirling, Z. Phys. C53 (1992) 629. 

[41] S. Catani, Yu.L. Dokshitzer, F. Fiorani and B.R. Webber, Nucl. Phys. B377 (1992) 
445. 

[42] L3 Collaboration, B. Adeva et a/., Phys. Lett. B248 (1990) 464; 

OPAL Collaboration, M.Z. Akrawy et al. , Phys. Lett. B235 (1990) 389; 

OPAL Collaboration, M.Z. Akrawy et al, Z. Phys. C49 (1991) 375; 

ALEPH Collaboration, D. Decamp et al. , Phys. Lett. B255 (1991) 623; 

DELPHI Collaboration, P. Abreu et al, Phys. Lett. B247 (1990) 167; 

DELPHI Collaboration, P. Abreu et al, Z. Phys. C54 (1992) 55; 

ALEPH Collaboration, D. Decamp et al. , Phys. Lett. B284 (1992) 163; 

OPAL Collaboration, P.D. Acton et al, Z. Phys. C55 (1992) 1; 

ALEPH Collaboration, D. Buskulic et al, Z. Phys. C55 (1992) 39; 

ALEPH Collaboration, D. Buskulic et al, Z. Phys. C55 (1992) 209; 

OPAL Collaboration, P.D. Acton et al, Z. Phys. C59 (1993) 1. 

[43] See for a review: 

Chapter 12.5 of Ref. [|l| (and references therein). 

[44] L3 Collaboration, B. Adeva et al., Phys. Lett. B263 (1991) 551; 

DELPHI Collaboration, P. Abreu et al., Phys. Lett. B274 (1992) 498; 

OPAL Collaboration, G. Alexander et al., Z. Phys. C52 (1991) 543. 

[45] L3 Collaboration, B. Adeva et al, Phys. Lett. B248 (1990) 227; 

OPAL Collaboration, M.Z. Akrawy et al, Z. Phys. C49 (1991) 49. 

[46] DELPHI Collaboration, P. Abreu et al., Phys. Lett. B255 (1991) 466; 

ALEPH Collaboration, D. Decamp et al, Phys. Lett. B284 (1992) 151; 

DELPHI Collaboration, P. Abreu et al, Z. Phys. C59 (1993) 357; 


L3 Collaboration, B. Adeva et al, Phys. Lett. B248 (1990) 227; 

OPAL Collaboration, R. Akers et al, Z. Phys. C65 (1995) 367; 

K. Hamacher, Proceedings of the XXVIIIth International Conference on High Energy 
Physics, Warsaw, Poland, 25-31 July 1996, page 749. 

[47] ALEPH Collaboration, R. Barate et al, Z. Phys. C76 (1997) 1. 

[48] L3 Collaboration, B. Adeva et al, Phys. Lett. B248 (1990) 227; 

OPAL Collaboration, M.Z. Akrawy et al, Z. Phys. C49 (1991) 49. 

[49] L. Lonnblad and T. Sjostrand, Phys. Lett. B351 (1995) 293. 

[50] G. Gustafson, U. Pettersson and P.M. Zerwas, Phys. Lett. B209 (1988) 90; 

T. Sjostrand and V.A. Khoze, Phys. Rev. Lett. 72 (1994) 28; Z. Phys. C62 (1994) 281; 
G. Gustafson and J. Hakkinen, Z. Phys. C64 (1994) 659. 

[51] Proceedings of the Workshop ‘Physics at LEP2’, eds. G. Altarelli, T. Sjostrand and F. 
Zwirner, CERN Report 96-01. 

[52] Proceedings of the ‘LEP2 Phenomenology Workshop’, Oxford, UK, 14-18 April 1997, 
J. Phys. G24 (1998) 271. 

[53] See, e.g.: 

Proceedings of the Workshop ‘e + e _ Collisions at 500 GeV. The Physics Potential’, Mu¬ 
nich, Annecy, Hamburg, ed. P.M. Zerwas, DESY 92-123 A/B/C, 1992-1993; 
Proceedings of the Workshop ‘Physics with e + e _ Colliders’, Annecy, Gran Sasso, Ham¬ 
burg, ed. P.M. Zerwas, DESY 97-100, May 1997. 

[54] T. Ohl, Proceedings of CRAD’96, Third International Symposium on Radiative Cor¬ 
rections, Cracow, 1-5 August 1996, Acta Phys. Pol B28 (1997) 847. 

[55] E. Accomando, A. Ballestrero and M. Pizzio, Nucl. Phys. B512 (1998) 19. 

[56] F. Yuasa and Y. Kurihara, in preparation. 

[57] G. Montagna, M. Moretti, O. Nicrosini and F. Piccinini, Eur. Phys. J. C2 (1998) 483. 

[58] E. Accomando, A. Ballestrero and E. Maina, Phys. Lett. B362 (1995) 141. 

[59] F.A. Berends, W.L. van Neerven and G.J. Burgers, Nucl. Phys. B297 (1988) 429; 
Erratum, ibidem B304 (1988) 95; 

E.A. Kuraev and V.S. Fadin, Sov. J. Nucl. Phys. 41 (1985) 466; 

G. Altarelli and G. Martinelli, Proceedings of the Workshop ‘ Physics at LEP', eds. J. 



Ellis and R. Peccei, Geneva, 1986, CERN 86-02; 

R. Kleiss, Nucl. Phys. B347 (1990) 29; 

O. Nicrosini and L. Trentadue, Phys. Lett. B196 (1987) 551; Z. Phys. C39 (1988) 479. 

[60] G. Montagna, O. Nicrosini and F. Piccinini, Phys. Lett. B406 (1997) 243. 

[61] See, e.g.: 

C. Itzykson and J.-B. Zuber, “Quantum Field Theory” (McGraw-Hill, 1985). 

[62] M.Jacob and J.C. Wick, Ann. Phys. 7 (1959) 404. 

[63] J.D. Bjorken and M.C. Chen, Phys. Rev. D154 (1966) 1335; 

O. Reading-Henry, Phys. Rev. 154 (1967) 1534. 

[64] P. De Causmaecker, R. Gastmans, W. Troosts and T. T. Wu, Phys. Lett. B105 (1981) 
215; Nucl. Phys. B206 (1982) 53; 

F.A. Berends, R. Kleiss, P. de Causmaecker, R. Gastmans, W. Troosts and T. T. Wu, 
Nucl. Phys. B206 (1982) 61; ibidem B239 (1984) 382; ibidem B239 (1984) 395; ibidem 
B264 (1986) 243; ibidem B264 (1986) 265. 

[65] M. Caffo and E. Remiddi, Helv. Phys. Acta. 55 (1982) 339. 

[66] G.R. Farrar and F. Neri, Phys. Lett. B130 (1983) 109. 

[67] G. Passarino, Phys. Rev. D28 (1983) 2867; Nucl. Phys. B237 (1984) 249. 

[68] F.A. Berends, P.H. Daverveldt and R. Kleiss Nucl. Phys. B253 (1985) 441. 

[69] R. Kleiss and W.J. Stirling, Nucl. Phys. B262 (1985) 235. 

[70] R. Kleiss and W.J. Stirling, Phys. Lett. B179 (1986) 159. 

[71] J. Gunion and Z. Kunszt, Phys. Lett. B161 (1985) 333. 

[72] K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B274 (1986) 1. 

[73] C. Mana and M. Martinez, Nucl. Phys. B287 (1987) 601. 

[74] A. Ballestrero and E. Maina, Phys. Lett. B350 (1995) 225. 

[75] Z. Xu, Da-Hua Zhang and L. Chang, Nucl. Phys. B291 (1987) 392. 

[76] H. Murayama, I. Watanabe and K. Hagiwara, HELAS: HELicity Amplitude Subrou¬ 
tines for Feynman Diagram Evaluations, KEK Report 91-11, January 1992. 



[77] S. Moretti, Phys. Rev. D50 (1994) 2016. 

[78] D. Zeppenfeld, Int. J. Mod. Phys. A3 (1980) 2175. 

[79] P. Cvitanovic, Phys. Rev. D14 (1976) 1536. 

[80] G.P. Lepage, Jour. Comp. Phys. 27 (1978) 192. 

[81] JADE Collaboration, W. Bartel et al., Phys. Lett. B123 (1983) 460; Z. Phys. 33 (1986) 
23. 

[82] Yu.L. Dokshitzer, contribution cited in Report of the Hard QCD Working Group, Proc. 
Workshop on Jet Studies at LEP and HERA, Durham, December 1990, J. Phys. G17 
(1991) 1537. 

[83] Yu.L. Dokshitzer, G.D. Leder, S. Moretti and B.R. Webber, J. High En. Phys. 8 (1997) 
1 . 

[84] G.D. Leder, Nucl. Phys. B497 (1997) 334. 

[85] S. Bethke, Z. Kunszt, D.E. Soper and W.J. Stirling, Nucl. Phys. B370 (1992) 310. 

[ 86 ] A. Signer and L. Dixon, Phys. Rev. Lett. 78 (1997) 811. 

[87] For a review, see, e.g.; 

S. Bethke, in Proceedings of the 42nd Scottish LIniversities Summer School in Physics 
on ‘High Energy Phenomenology’, August 1-21 1993, St. Andrews, UK. 

[ 88 ] T. Sjostrand, Coinp. Phys. Commun. 39 (1984) 347; 

M. Bengtsson and T. Sjostrand, Comp. Phys. Commun. 43 (1987) 367. 

[89] G. Marchesini, B.R. Webber, G. Abbiendi, I.G. Knowles, M.H. Seymour and L. Stanco, 
Comp. Phys. Commun. 67 (1992) 465. 

[90] T. Sjostrand, Comp. Phys. Commun. 82 (1994) 74. 

[91] R. Odorico, Comp. Phys. Commun. 59 (1990) 527. 

[92] L. Lonnblad, Comp. Phys. Commun. 71 (1992) 15. 

[93] G. Cowan, in Ref. |72| . 

[94] S. Moretti and W.J. Stirling, preprint RAL-TR-1998-037, DTP/98/32, Cavendish- 
HEP-97/12, July 1998. 


[95] J. Andre and T. Sjostrand, Phys. Rev. D57 (1998) 5767; 

J. Andre, preprint LU-TP-97-12, June 1997, hep-ph/9706325. 

[96] G. Marchesini, B.R. Webber, I.G. Knowles, M.H. Seymour, G. Corcella, S. Moretti and 

K. Odagiri, in preparation. 

[97] E. Fahri, Phys. Rev. Lett. 39 (1977) 1587. 

[98] Mark-J Collaboration, D.P. Barber et al, Phys. Rev. Lett. 43 (1979) 830. 

[99] G. Parisi, Phys. Lett. B74 (1978) 65; 

J.F. Donoghue, F.E. Low and S.Y. Pi, Phys. Rev. D20 (1979) 2759. 

[100] T. Chandrahoman and L. Clavelli, Nucl. Phys. B184 (1981) 365. 

[101] J.D. Bjorken and S.J. Brodsky, Phys. Rev. D1 (1970) 1416. 

[102] S. Catani and B.R. Webber, J. High En. Phys. 10 (1997) 5; Phys. Lett. B427 (1998) 
377. 


[103] For a review, see |K| and references therein. 


[104] ALEPH Collaboration, ‘QCD studies with e + e~ annihilation data at 130, 136, 161 
and 172 GeV’, contribution to EPS-HEP97, Jerusalem, 19-26 August 1997, reference 
number 629. 


[105] G. Gustafson, U. Pettersson and P.M. Zerwas, Phys. Lett. B209 (1988) 90; 

T. Sjostrand and V.A. Khoze, Phys. Rev. Lett 72 (1994) 28; Z. Phys. C62 (1994) 281; 
G. Gustafson and J. Hakkinen, Z. Phys. C64 (1994) 659; 

L. Lonnblad, Z. Phys. C70 (1996) 107; 

J. Ellis and K. Geiger, Phys. Lett. B404 (1997) 230. 

[106] B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rep. 97 (1983) 
31. 


[107] T. Sjostrand and V.A. Khoze, in Ref. [ pR5 


[108] ‘Determination of the Mass of the W Boson’, Z. Kunszt et a/., in ‘Physics at LEP2’, 
eds. G. Altarclli, T. Sjostrand and F. Zwirner, CERN Report 96-01, Vol.l, p.141-205 
(1996). 


[109] See, e.g.: 

B.R. Webber, in Ref. 


52 

















[110] P. Igo-Kemenes (convener) and J.H. Kiihn (advisor), in Ref. [Q, part A, page 255. 

[111] P. Igo-Kemenes (convener) and J.H. Kiihn (advisor), in Ref. |nj, part A, page 327. 


[112] T.L. Barklow, P. Bussey, P. Chen, M. Drees, R.M. Godbole, W. Kozanecki, F. Richard 
and J.C. Sens, in Ref. |)3|, part B, page 837. 

[113] T. Barklow, P. Chen and W. Kozanecki, in Ref. [Q, part B, page 845. 

[114] O. Nicrosini and L. Trentadue, in Ref. [^9 . 


[115] See, e.g.: 

Proceedings of the Workshop on Photon Radiation from Quarks, S. Cartwright ed., 
Annecy, France, 2-3 Dec. 1991. CERN 92-04; 

OPAL Collaboration, G. Alexander et a/., Phys. Lett. B264 (1991) 219; 

OPAL Collaboration, P.D. Acton et al, Z. Phys. C 54 (1992) 193; 

ALEPH Collaboration, D. Decamp et al., Phys. Lett. B264 (1991) 476; 

DELPHI Collaboration, P. Abreu et al, Z. Phys. C 53 (1992) 555; 

L3 Collaboration, O. Adriani et al ., Phys. Lett. B292 (1992) 472. 


[116] E.W.N. Glover and A.G. Morgan, Z. Phys. C62 (1994) 311; 

A. Gerhmann-De Ridder and E.W.N. Glover, Nucl. Phys. B517 (1998) 269. 


[117] J. Hakkinen and H. Kharraziha, Comput. Phys. Commun. 100 (1997) 311. 




















6 Appendix I 


In this Appendix we present the analytical expressions for the hclicity amplitudes relevant to 
the three processes (,!])— (;3|). As already mentioned, we make use of the formalism introduced 
through Refs. |j5], |73], Each amplitude will be expressed in terms of a suitable 

combination of the X and Z functions defined in Subsect. |2T] (and of scalar products 
involving the external four-momenta). The directions of the latter are such that the both 
the initial state electrons and positrons and all the partons q, q', q", q, q', q"g in the final 
states are considered outgoing. 

In order to minimise the length of the formulae, we adopt the following short-hand 
notations. The index (a) (i.e., in round brackets) stands for the momentum p a and the 
corresponding polarization A a of any real particle a — 1, ... 8 . In the case where external 
gluons are involved, additional indices will appear, such as [a] (i.e., in squared brackets), 
to which are associated the auxiliary momentum q a of the gluon a and its hclicity A a (see 
eq. (|T6|)). The notation {k} (i.e., in curly brackets) refers to both the momentum p and 
the internal hclicity A associated to a propagator. We further make the convention that 
all indices {k} are summed over, and that the sum extends to all possible values of all 
internal hclicities, A, A',... = ±, as well as to the momenta appearing in the last lines of the 
forthcoming expressions for the amplitudes of processes ( 0 )^( 0 ) (in the spirit of eq. ( 0 ))Q. 

To each of the basic amplitudes are associated propagators functions. In the case of 
off-shell fermions, they have the form 

1 


Dm) = 


(E iPi 


m 


2 ’ 
/ 


(30) 


whereas in the case of bosons one gets 


Dv({l}) (Z,Pi) 2 -M*+iM v r v 


(31) 


In eqs. (|30|)-(|3T|), / is the flavour g, q' or q" of a virtual fermion line, whereas V = 7 , Z or 
g, being My = Ty = 0 if V = 7 or g. In order to assign the correct fermionic propagator 
functions to the various amplitudes, one will have to multiply the latter for as many terms 
of the form ( |30| ) as the number of indices {i},... appearing in the last line of the amplitude 
equations, with the sums extending to all the values assumed by i. A subscript of the 
form {...}/ will indicate the flavour f = q, qq" to be assigned to the off-shell fermion line. 
As for the bosonic propagators in the case V = g (i.e., for gluons), these will only appear 
in the functions A, B or C (see below). Their appropriate expressions will be given while 
defining the latter. Finally, as the incoming V = 7 , Z currents are a common feature to all 


16 Note the order of the indices in the sequence {a, b 1 c, ...} has no importance. 










diagrams, we will factorise them in the expressions of the amplitudes, an operation which 
will be recalled by the use of the subscript V. Given these simple rules, we believe it trivial 
to derive the correct expressions for all propagators, so these will not appear explicitly in our 
formulae for the amplitudes. In the latter, we will present the expressions for the topologies 
of Fig. [|, in the form Ty(Jl, {A}), with i — 1, ...8 for process (P, i — 1, ...10 for process (||), 
and i = 1, ...3 for process (|3|), and where II represents (a product of) permutations n of 
gluon colour indices and {A} refers to the hclicities of the external particles. 

As for the couplings, the notation fV refers to the pair of chiral indices (cr, Cl) entering 
in the expressions given in Tabs. M and associated with the vertex involving a fermion 
/ = e,q,q',q" and a gauge vector V = 7 ,Z,g, according to Tab. f|. In particular, the 
notation qg is meant to describe both a true and a fictitious quark-gluon vertex, the latter 
being used to define the polarization of the gauge particle (for which (cr, Cl) = (1,1)), 
according to eq. m- 

Graphically, a Z function is represented as (hereafter, one symbolically has that () = (), [] 
or {}) 


Z (( a ); (b ); (c); (d ); fV ; f'V) = 



(32) 


On the same footing as was done in Ref. 


a 


we found it convenient to introduce 


the function describing a three-gluon vertex, here generalised to the case of off-shell gauge 


vectors, each connected to a fermion limp], 


17 Here and in the following, we write the momentum of the gluon attached to the fermion line labelled 
by ab in terms of those of the others involved. This convention will be helpful both in writing the Feynman 
diagrams in a more compact analytic form as well as in calculating them numerically in a more efficient 
manner, as this way the labels a and b only appear as arguments of the A, B and C functions, and not in 
the internal momentum summations. In fact, the line ab will eventually be identified with the off-shell one 
connected to the incoming 7 , Z current in Fig. |]ri-c. 







= Z{(a)\{b)\{c)\(d)\qg\qg) 
x 2X((e ); {z}; {f);qg) 

- Z((c)-,(d)](e)-,(f)-,qg;qg) 
x X((a)]{i};(b);qg) 

— same[(c, d) <->• (e,/)]. (33) 


To each A term is always associated the propagator D g ({i} + {j}) = 1 /{J2iPi + HjPj) 2 - 
Note that the summation {i} extends to % — c, d if (d) 7 ^ [d], or to i — c otherwise and 
similarly for {)} with respect to the indices e,/. In addition, if (d) 7 ^ [d]((f) 7 ^ [/]), the 
additional term D g ({i}){D g ({j})) appears. 

In addition to eq. (|33]), we also need the expression describing two connected three-gluon 
vertices and the associated four-gluon vertex with same colour structure. Assuming again 
that the gluon propagator connected to ab is expressed in terms of the momenta associated 
to the legs c, one has 


B((a ); ( b ); (c); (d ); (e); (/); (g ); (h)) 



c b c b 



= Z((a) ; (b) ; (c); (d) ; gg; qg) 

x [ 2 ( x ((g ); (0; ( h ); gg) + A "((^); 0'}; ( h ); <?#)) x 

2 X((e);{i};(/);g^)] 

+ Z((a) ; ( 6 ); (g); (/i); gg; qg) Z((c) ; (d); (e); (/); gg; qg) x 

( 2 5>fc + 53+ HPj) ' CE 'Pi) 

k i j i 

+ Z((a) ; ( 6 ); (e); (/); gg; gg) Z((c); (d); (g); (/i); gg; qg) x 
* j 

~ ^((c); (d); (e); (/) ;gg;gg) 

x ( 2 (^((g); (0; W; <?#) + ^((#); 0'}; W; gg)) x ((a ); {*}; ( & ); Qg) + 

(X((a) ; {k}; (b) ; gg) - X((a ); {j}; (6) ; gg) - A"((a); {i}; (b) ; qg))X((g ); {i}; (/i); gg)] 
+ 2X((e);{i};(/);gg) 
x [-Z((a);(b);(g);(h);qg;qg) 

x ( 2 X((c); {k}; ( d ); gg) + X((c); {j}; (d); gg)) 

+ ^((c); <d); (g); (h) ;gg;gg) 

x (A"((a); {/c}; ( 6 ); gg) - X((a) ; {j}; ( 6 ); gg) - -A((a); {i}; (b) ; gg))] 

— same[(c, d) <->• (e, /)] (34) 


The propagators associated to 5 functions assume the form + {j} + {/c})i4 s ({i} + 



0'}) = 1 /(EiPi + J2jPj + E kPk) 2 /{J2iPi + T,jPj) 2 - The summations {*}, {j} and {k} are 
over the indices c,d, e, f and g,h if ( d) ^ [d], (/) ^ [/] and (h) ^ [h], or c, e and g 
otherwise, respectively. Moreover, as previously explained, if (d) ^ [d]((f) ^ [f])[{h) ^ [/i]], 
the additional term D g ({i})(D g ({j}))[D g ({k})] has to be inserted. 

Finally, we further need the function describing three connected three-gluon vertices and 
the two associated graphs with one four-gluon and one three-gluon vertex with identical 
colour matrices, 




= {(2 X((i) ; {k}] (j) ; gg) + 2 X{(i) ; {/}; {j ); qg ) + 2X((i) ; {m}; (j) ; gg) + X((i) ; {ra}; (j) ; qg)) 

x[ 2 (x((g) ; W; W ;gg) + x ((g ); (0; W ;gg)) x 
(2X((e); {/c}; (, f);qg)Z((a ); (&); (c); (d) ;gg;gg) - 
X({a ); {&}; ( 6 ); gg) Z((c); (d); (e); (/); gg; gg)) 

+ ( A "((°); M; (6); gg) - ^((a); W; ( & ); gg) - A "((°); W; (&); gg)) x 

( 2 X((e); {/c}; (/); gg) Z((c); (d); (g); (h ); gg; gg) - 
^((g>; W; ( h ); gg) ^((c); (d); (e); (/); qg-, qg)) 

- ( 2X ((c ,); M; (d); gg) + ^((c); {i}; <d); gg)) x 

2X((e); {A;}; (/); gg) Z((a); (6); (g); (/i); gg; gg) 

+ (2Zp™ + Zp* + J2Pk) ■ (Zp*) x 

m l k k 

Z{{a ); (b ); (g) ; ( h ); gg; gg) Z((c); (d); (e); (/); gg; gg)] 

- same[(a, 6) <-»• (i, j); {k} + {/} + {m} + {n} <-> — {n}]} 

+ ^((«); (&); (*); O’) 5 gg; gg) Z( c ); (d>; (e); (/); gg; gg) 

x[2(^«g) ; {&}; (h);qg)+X((g ); {/}; (h);qg)) x 

(ZPfc + ZP* + ZPm + 2 ZPn) • (Zp*) 

fe l m n k 

+ ^((g>; W; W;gg) x 

(Zp* + Zp* + Zp™ + 2 Zp«) ■ (Zp™ - J2pi - Zp*) 

k l m n m l k 

~ ( x ({g ); {k}] (h ); gg) + X((g) ; {/}; (h ); gg) + 2X((g ); {n}; (/i); gg)) x 

(Z Pfc + Z Pi + 2 Z• (ZPfc)] 

k l rn k 

+ 2 Z((a); ( 6 ); (i ); (j); gg; gg) X((e); {&}; (/); gg) 
x[Z«c); (d); (g); (./i) ;gg;gg) 

(ZPfc + Zp' + ZPm + 2 Zp«) ‘ (ZPfc + ZPi.^ZP™) + 

k l rn n k l rn 

( 2 A"((c); {m}; (d); gg) + X((c); {/}; (d); gg)) x 

( A "((g) ; W; (M ; gg) + x ((g) ; W; (M; gg) + 2X((g) ; {n}; (h ); gg)) — 

2(2X((c); {n}; (d); gg) + X((c); {m}; (d); gg) + A"((c); {/}; (d); gg))) x 

(Z(g); W; W ;gg) + ^((g); W; W; gg))] 

+ (Zp^ + Zw + ZPm) 2 x 

k l rn 

{Z((a ); ( b ); (g); (h);qg;qg) x 
[2Z«c) ; (d) ; (i); (j); gg; gg) A"((e) ; {k}; (/); gg) - 



Z((c) ; (d ); (e) ; (/) ; gg; qg) X((i) ; {k}; (j) ; gg)] 

-Z({g ); </i); (?); (/) ;gg;gg ) x 
[2Z((a); ( 6 ); (c); (d); gg; qg) X((e ); {fc}; (/); qg) - 
Z«c) ; (d ); (e) ; (/) ; qg] qg) X((a) ; {k}] ( b ) ; gg)]} 

+ (X)p* + X)pj) 2 ^(( c );( d >;(s>; W;gg;gs) x 

k l 

[2Z((a) ; ( 6 ); (e); (/); gg; qg) x 

(X((i); {A;}; (j); gg) + X((i); {/}; (j ); gg) + X((i ); {m}; (j); qg)) + 

z (( e ); (/); (0; {j)\gg\gg) x 

(X((a); {n}; (6); qg) - X((a) ; {m}; (6); gg) - 
A "((«); {^}; ( b ); gs) - ^((a); {k}] (b ); gg)) - 
Z((a); ( 6 ); (i ); (j); gg; gg) x 

(2X((e); {n}; (/); gg) + X((e); {m}; (/); gg) + X((e);{/};(/); gg))] 

— same[(c, d) (e, /)]. (35) 

The propagator functions associated to these terms are of the form D g ({k} + {/} + {m} + 
{n})D g {{k} + {1} + {m})D g ({k} + {/}) = 1/(E kPk + E 1 P 1 + E m Pm + E n Pn) 2 /(EkPk + 
EiPi + EmPm) 2 /(EfcPfc + EiP«) 2 - The summations {k}, {/}, {m} and {n} are here over c,d, 
ej, g, h and i,j if (d) ^ [d\, (/) ^ [/], (h) j- [h] and (j) ^ [j], or c, e, g and i otherwise, 
in the same order. Also in this case, terms of the forms D g ({k}), D g ({l}), D g ({m}) and 
D g ({n}) will appear if (d) ^ [d], (/) ^ [/], (h) ^ [h\ and/or (j) \j], respectively. 

Note that further manipulations would be possible, though in the form given in eqs. 


(134| ) and (|35|) , the functions A, B and C are suitable for immediate implementation in a 
numerical program, as we have done in our calculations. The expressions (|32|), (p3|), (|34|) 
and (^5|) are now the basic ingredients needed to describe the spinor part of all the topologies 
appearing in Fig. [|a c. 

In the remainder of this Appendix, we present the hclicity amplitudes associated to the 
typical topologies in Fig. [|a-c, stripped off their colour structure (operation that we recall 
by means of the introduction of the argument II). For sake of illustration, we assume that 
the labels a and b always refer to the quarks and antiquarks of flavour q in all reactions 
(0M0), that c and d refer to gluons in reaction (|Tj) and to the quarks and antiquarks of 
flavour q' in reactions (J2])-([|) and that e and / refer to gluons in reactions (|I])-(@) and 
to the quarks and antiquarks of flavour q" in reaction (|3]). The permutations chosen as 
example for the gluons labelled cdef in reaction ([[]) and ef in reaction (| 2 |) as well those 
for the internal coloured gauge vectors in processes (|2|)-(|3|) ought to be self-evident in the 
forthcoming helicity amplitudes. 





T i y ( n >( A }) = Z{{a)\{i}\{c)\[(^\qg\qg) 
x Z({i};{j}-,(d)’,[d\;qg;qg) 
x Z{{j}-{k} ] {l)-{2)-qV-eV) 
x Z({h};{l};(e);[e]\qg;qg) 
x Z{{l}\{b)\{f)\[f]\qg\qg), 

(0 = {a, c} g , {j} = {a, c, d } q , {/c} = { 1 , 2 , a, c, rf} q , 

{/} = {1,2, a, c, d, e} q . (36) 

T 2 y ( n , {A}) = 4((a);{i};(c);[c];(d);[d]) 
x Z({i}-{j}-(l)-(2)-qV-eV) 
x Z({j}]{k}](e)-,[e\,qg\qg) 

x ^(W;(6);(/);[/];gp;g^), 

{*} = {a, c, d} g , {j} = {1, 2, a, c, d}„ {k} = {1, 2, a, c, d, e} q . (37) 

T 3 V (H, {A}) + Tf (n, {A}) = B((a);{z};(c);[c];(d);[d];(e);[e]) 

x Z({z};{j};(l);(2);qV;eV) 

x Z ({J}; (&); (/);[/]; w, qg), 

{*} = {a, c, d, e} q , {j} = {1, 2, a, c, d, e} q . (38) 

T 5 (H, {A}) = ^((a);{i};(c);[c];(d);[rf]) 
x Z({i};{j};(l);( 2 );^;ey) 
x A( 0 }; (6); (e); [e]; (/);[/]), 

{*} = {a, c, {j} = {1, 2, a, c, d} q . (39) 

Tq (n, {A}) + (n, {a}) + r 8 y (n, {A}) = r((a) ; {<} ; ( i)-(2)- q v-eV) 

X C'd*}; ( 6 ); (c); [c]; (d); [d]; (e); [e]; (/); [/]), 
{*} = {!, 2, a},. (40) 





T i (n, {A}) = Z{{a)\{i}\ (c); (d);qg;qg) 
x Z({t}-{j}-{l)-(2)-qV-eV) 
x A({j}; (6); (e);[e] ; (/);[/]), 

{*} = {a, c, d} g , {j} = {1, 2, a, c, d} g . (41) 

r 2 V (n, {A}) = Z((a); {i}; (c); (d);qg;qg) 
x Z({i};{j};(l);(2);gy;ey) 

x Z({k};(b)-,(f);[f]-,qg-qg), 

{*} = {a, c, d} g , {j} = {1, 2, a, c, d} 9 . {j} = {1, 2, a, c, d, e} q . (42) 

T 3 V (n,{A}) = Z((a)]{i};{k};(d);qg-,qg) 
x Z((c);{/c};(e);[e];g^;g^) 
x Z({i};{i};(l);(2);gy;el/) 
x Z({j}-,(by,(f);[f]-,qg-,qg), 

{*} = {a, c, d, e} q , {j} = {1, 2, a, c, d, e} q , {/c} = {c, e},/. (43) 

T/(n,{A}) = A((a);{i};(c);(d);( e );[ e ]) 
x Z({i};{j};(l);(2);gy;ey) 

x ^(0'};(6);(/);[/];e^;?^), 

{*} = {a, c, d, e} q , {.)} = {1, 2, a, c, d, e} g . (44) 

T^(n,{A}) = Z((a);{i};(l);(2); g V;eV) 
x Z({i};(b);{k}](d)]qg;qg) 
x Z((c); {j}; (e); [e];qg; qg) 

x £(0'};{fc}; (/); [/];?s;gs), 

{*} = {1, 2, a } q , {j} = {c, e}g/, {&} = {c, e, f} q >. 


( 45 ) 





T e ( n , {A}) = Z((a);{t};(l);(2); g V;eV) 
x Z{{i}](b)]{j}\{d)\qg\qg) 
x ^((c);{j};(e);[e] ; (/);[/]), 

{*} = {1, 2, a} q , {j} = {c, e, /}y. (46) 


T 7 (n, {A}) + r 8 y (n, {A}) + T 9 y (n, {A}) = Z((a);{i};(l);(2); g V; e V) 

x Sdi}; (6); (c); (d); (e); [e]; (/); [/]), 

{*} = {!, 2, a},. (47) 


T io(n,{A}) = Z((a);{i};(l);(2); g V;eV) 

X A({»}; (6); (j};(d); (/);[/]) 

X ^((c);0'};(e);[e];g^;g^), 

{*} = {1, 2, a} g , (. 7 } = {c, e} g /. (48) 


e + (pi, Ai) + e (p 2 , A 2 ) -> g(p a , A a ) + g(p b , A 6 ) + g'(p c , A c ) + q'(p d , X d ) + q"(p e , A e ) + q"(p f ,X f ) 

T i (n, {A}) = Z((a); {i}; (1); (2); qV ; eV) 
x Z({i};(b);{j}](d)-,qg-,qg) 
x Z((c)]{j};(ey,(f);qg;qg), 

{*} = {1, 2, a} q , {j} = {c, e, f} q >. (49) 


r 2 V ( n , {A}) = Z((a);{i}](c);(d)- 1 qg-,qg) 
x Z{{i} ] {j}-{l)-{2)-qV-eV) 
x Z({j};(b)](e);(f);qg-,qg), 

{*} = {a, c, d } q , {j} = {1, 2, a, c, d} q . 


( 50 ) 






^CMA}) = Z((a);{i};(l);( 2 );gF;eF) 
x A ({ih(b)](c);(d);(e);(f)), 

{i} = {1,2 ,a} q . (51) 


7 Appendix II 


In this Appendix we describe how we have proceeded to obtain the QCD colour factors^ 
relevant to the three reactions (HD-©. In what follows, we assume that all internal and 
external colours are always summed over. 

The amplitudes squared for processes ©)-©, summed/averaged over the final/initial 
colours and spins, can be written in the general form 




n N <) (n —) E E V(n, {A»M'(n', {A})c n - n \ 

{A} n,n' 


n p \ 


(52) 


The factor 1/4 accounts for the average over the initial spins whereas the EM and QCD 
couplings are obtained via the relations e = \/Aira em and g s = >/ 47 tck s (in natural units) from 
the constants introduced in Subsect. |2.3| . In the above equation IT^ represents (products 
of) permutations of the gluon indices, {A} indicates the helicities of all external particles 
Aj (with i = 1, ... 8 ), M(II('),{A}) is the coefficient of the colour matrix Tj 1<! (where / 
represents a sequence of flavours ij...), that is, the full spinor amplitude (here, {i} = { 1 , 2 }) 


M(n«, {A}) 


E E ( E n-in't {A})£V(W) 

Tt(qq'(q")) i= 1 \V =7 ,Z 


(53) 


with n — 8 , 10 and 3 for process ([!]), (Q) and ©, respectively, and is the appropriate 

colour factor, C n,n ' = TrlTj^T^ Q. Here, Z^Tr^qg")) refers to the sum over all possible flavour 
permutations in events with at least four quarks. Furthermore, n p \ is a statistical factor 
corresponding to each ?r p -tuple of identical final state particles p appearing in reactions 
(IB whereas N t is the normalisation of the fermion current associated to the gluon i, see 
eq. (pT|). In reaction © P = g and n g = 4, for © one has p = q,q,g with n q = n q = 1(2) 
if q 7 ^ q\q = q') and n g = 2, whereas for © one gets p = q,q with n q = n q = 1(2) [3] if 
g / / / q'\q = q' ^ q")[q = q' = q"}. As for the gluon normalisations, the product ]ffi 
extends to i — 5, 6 , 7, 8 in (|I|), to i — 7, 8 in ©) whereas it does not appear for ©• 


18 In order to check their correctness we have calculated the traces of the QCD matrices both using an 
analytical procedure as well as the computer program COLOUR [117]. 


19 Therefore, our multidimensional integrations must cover the whole six-particle phase space. 









In order to obtain the colour factors C n,n ' we have proceeded as follows. First, we have 
decomposed a four-gluon vertex into three three-gluon vertices, each with a different colour 
and kinematic structure in factorised form. One such decompositions, though not unique, 
that we have adopted here, is the one pictured in Fig. [L7| (from Ref. ||117| ). Then, we have 
repeatedly exploited the equation 


A ,B-j _ .A.B +B ±A • rABC + C 

' i " \ij — ^ik^kj ^ik^kj ‘■J "ij 


(54) 


relating the matrices tf- (with A = 1,...8 gluon and i,j = 1, ...3 quark colours) of the 
fundamental representation to the structure constants f ABC of the SU(3) gauge group^], in 
order to rewrite the colour parts of diagrams involving triple and quadruple gluon vertices 
as a combination of the basic colour structures typical of the ‘abelian graphs’, i.e., those 
non involving gluon self-couplings. If one does so, than the colour piece of each diagram can 
be expressed as an appropriate combination of different permutations of products involving 
tfj matrices only^J 

In case of process ([!]), all colour structures in Fig. 0 a are linear combinations of terms 
of the form (here, II = tt and / = ij) 


mmsmg) = (^W)* 8678 ’, 


(55) 


where A,B,C,D(i,j) are the gluon(quark) colour and 7 t( 5678) is one the possible permu¬ 
tations of the gluon indices A, ...D = 5, ... 8 . For four external gluons, one has 4! of these, 
that is 24, twelve of which will be cyclical and twelve anticyclical, so that the number of 
fundamental colour factors C ,7r ’ 7r, needed to describe the colour content of process ([[]) is 576 
in total. These are given in Tab. || for the cyclical 7 r'-permutations only, as the remaining 
colour factors can be obtained by symmetry. 

In case of process (2J), the colour structures appearing in Fig. |3|b can be rewritten as 
linear combinations of terms of the four forms (here, II = in the hrst(last) two cases 

and / = ijkl ) 


Ti ijkl (qqq'q'gg) = (l¥l B )f 78) (l X ) w , 

(56) 

TZ ijk i(qqq'q'gg) = (t x )ij(t x t A t B )$ xn) 

(57) 

Tr’^kMqq'q'gg) = (t x t% {X7 \t x t B )l? X8 \ 

(58) 


20 As a possible representation of the generators we have used the Gell-Mann one. 

21 Recall that if more than one coloured fermion line is present, such a combination is not unique, though 
the final answer obviously is. 




Tr [7% ( gggggg) (T%' (gggggg )) f ] 


£^7T,7T ' 

5678 

5786 

5867 

6758 

6875 

6587 

7856 

7568 

7685 

8576 

8657 

8765 

5678 

256 

4 

4 

4 

-5 

4 

71/2 

4 

-5 

-5 

-5 

-14 

5786 

4 

256 

4 

-5 

-14 

-5 

4 

4 

-5 

4 

71/2 

-5 

5867 

4 

4 

256 

71/2 

-5 

4 

-5 

-5 

-14 

4 

4 

-5 

6758 

4 

-5 

71/2 

256 

4 

4 

-5 

4 

4 

-14 

-5 

-5 

6875 

-5 

-14 

-5 

4 

256 

4 

-5 

71/2 

4 

-5 

4 

4 

6587 

4 

-5 

4 

4 

4 

256 

-14 

-5 

-5 

-5 

4 

71/2 

7856 

71/2 

4 

-5 

-5 

-5 

-14 

256 

4 

4 

4 

-5 

4 

7568 

4 

4 

-5 

4 

71/2 

-5 

4 

256 

4 

-5 

-14 

-5 

7685 

-5 

-5 

-14 

4 

4 

-5 

4 

4 

256 

71/2 

-5 

4 

8576 

-5 

4 

4 

-14 

-5 

-5 

4 

-5 

71/2 

256 

4 

4 

8657 

-5 

71/2 

4 

-5 

4 

4 

-5 

-14 

-5 

4 

256 

4 

8765 

-14 

-5 

-5 

-5 

4 

71/2 

4 

-5 

4 

4 

4 

256 

5687 

-32 

40 

-32 

-1/2 

-1/2 

-32 

31 

-1/2 

31 

-1/2 

40 

31 

5768 

-32 

-32 

40 

40 

31 

-1/2 

-1/2 

-32 

-1/2 

-1/2 

31 

31 

5876 

40 

-32 

-32 

31 

31 

-1/2 

40 

-1/2 

31 

-32 

-1/2 

-1/2 

6785 

-1/2 

31 

31 

-32 

-32 

40 

-1/2 

-1/2 

-32 

31 

-1/2 

40 

6857 

-1/2 

31 

40 

40 

-32 

-32 

31 

31 

-i/2 

-1/2 

-32 

-1/2 

6578 

-32 

-1/2 

-1/2 

-32 

40 

-32 

31 

40 

-1/2 

31 

-1/2 

31 

7865 

31 

-1/2 

31 

-1/2 

40 

31 

-32 

40 

-32 

-1/2 

-1/2 

-32 

7586 

-1/2 

-32 

-1/2 

-1/2 

31 

31 

-32 

-32 

40 

40 

31 

-1/2 

7658 

40 

-1/2 

31 

-32 

-1/2 

-1/2 

40 

-32 

-32 

31 

31 

-1/2 

8567 

-1/2 

-1/2 

-32 

31 

-1/2 

40 

-1/2 

31 

31 

-32 

-32 

40 

8675 

31 

31 

-1/2 

-1/2 

-32 

-1/2 

-1/2 

31 

40 

40 

-32 

-32 

8756 

31 

40 

-1/2 

31 

-1/2 

31 

-32 

-1/2 

-1/2 

-32 

40 

-32 

Table 6 

: The 

colour 

'actors 

£T7T,7r / 

multiplied by 

27, for 

e + e~ 

—*► qqgggg . Permutations n 


are along columns, permutations it' are along rows. Only half of the 24 V-permutations 
are given (cyclic ones), as the remaining colour factors can be obtained by symmetry: the 
(1 — 12) x (1 — 12) submatrix is identical to the (13 — 24) x (13 — 24) one (transposed), 
whereas the (13 — 24) x (1 — 12) submatrix is identical to the (1 — 12) x (13 — 24) one 
transposed. 













































































































































































































































































































































where A, B, X(i, j, k,l) are the gluon(quark) colour indices and 7d , )(A78), and 

(778) are again gluon index permutations. The symbol A" refers to any internal gluon 
A = 1, ...8 connecting the two quark-antiquark lines of colour ij and kl (and flavours q 
and q ', respectively). For two external gluons, one has 3! x 1! permutations of the former 
two kinds and 2! x 2! of the latter, for a total of 400 basic colour factors with 

m,n = 1,2, 3,4. We present these in Tabs. 0U for the combinations mn = 11,12,13, 33 
and 34. The colour factors ii\ 


are identical to the CA 


ones and the 


'' 22 ( 44 ) uu LHC V 11 ( 33 ) 

corresponding tables coincide with those numbered |7] and 0 here, respectively (tacitly 
assuming that the same gluon permutations are now made along different quark colour lines). 
Further notice that the colour factors C 2 3 can be obtained from Tab. |9] by exchanging 
7 with 8 in the tt' and tt" permutations, as can easily be understood by comparing the 
corresponding colour diagrams. Finally, the cases C^ 77 77 are identical to the C 2 ^ 77 ones 
whereas the C 2 ± 77 ’s coincide with the 17 ’s, so that the combinations mn = 14 and 24 
can both be deduced from Tab. |j. after appropriately assigning the quark colours ij and kl 
to the gluon permutations tt, tt’ and tt". 

Decomposing the colour structure of the Feynman diagrams in Fig. [3]b in terms of fl56[) - 
(|59|) is very convenient when dealing with diagrams which only differ in the exchange of the 
flavours, q q'. In fact, in this case, in order to obtain the colourless part of new diagrams 
one flips between themselves both the quark and antiquark labels in the expressions of the 
amplitudes (El -(^): i.e., a *-*■ c and b d (this way, both the momenta and hclicities 
get interchanged). As for the colour factors, these can easily be obtained from Tabs. 0-0 
if one notices that the colour matrices in ( 0 i are identical to those of eq. (j5B|) when one 
interchanges the (anti)quark colour labels, i.e., i <->■ k and j l, and so is the case for those 
in (|59| ) with respect to eq. (|58|). 

In case of identical flavours in process (Q), q — q', one further needs to calculate the 
diagrams which differ (in the spinor part) from the original ones of the case q 7 ^ q' in the 
exchange of the four-momentum and hclicity of one of the, say, quarks with those of the 
other, a c. Because of Fermi statistics, a minus sign factorises too. However, (anti)quark 
colour quantum numbers have to be exchanged as well. This introduces new colour struc¬ 
tures, identical to those in eqs. ([56|) - (|59|) except for the exchange i <-»■ k. Clearly, the colour 
factors building up the squared amplitude of the new set of graphs are identical to those 
already given in Tabs. M- In contrast, new ones are needed for the interference terms. 
These appear in the above tables in round brackets (note their different normalisation). 

In case of process (|3j), we recognise the following three basic colour structures, as the 


Tr [Tim ( qqq'q'gg ) { T u jk i ( qqq'q'gg )) f ] 

/''T7T,7r / 

0 11 

X78 

78X 

8X7 

87X 

X87 

7X8 

X78 

32 (-32) 

l/2(-l/2) 

l/2(-l/2) 

5(-5) 

-4(4) 

-4(4) 

78X 

l/2(-l/2) 

32(-32) 

1/2C-1/2) 

-4(4) 

5(-5) 

-4(4) 

8X7 

l/2(-l/2) 

l/2(-l/2) 

32(-32) 

-4(4) 

-4(4) 

5(-5) 

87X 

5(-5) 

-4(4) 

-4(4) 

32 (-32) 

l/2(-l/2) 

l/2(-l/2) 

X87 

-4(4) 

5(-5) 

-4(4) 

l/2(-l/2) 

32(-32) 

l/2(-l/2) 

7X8 

-4(4) 

-4(4) 

5(-5) 

1/2C-1/2) 

l/2(-l/2) 

32(-32) 


Table 7: The colour factors , multiplied by 9, for e + e _ —> qqq'q'gg (q ^ q'). In round 
brackets are given the additional ones needed in the case of identical flavours q — q' (for 
which the exchange of indices i <-> k is understood in one of the two interfering colour 
structures), multiplied by 27. Permutations i r are along columns, permutations n' are along 
rows. The first three are cyclical permutations, whereas the latter three are anticyclical. 
Note the symmetry in the colour factors: the (1 — 3) x (1 — 3) submatrix is identical to the 
(4 — 6 ) x (4 — 6 ) one (transposed), whereas the (1 — 3) x (4 — 6 ) submatrix is identical to 
the (4 — 6 ) x (1 — 3) one (transposed). 


fundamental ones to calculate all the colour factors pertinent to the graphs in Fig. |3]c (here, 
II = 7T and / = ijklmn ): 

T ltffcimn(9??WV') = (t X t Y )ij XY \t X )kl{t Y )mn, (60) 

T2- jk i rnn {qqq'q'q"q ") = {t X )ij(t x t Y )$ XY) (t y ) mn , (61) 

T £ j k l mn{qqq'q l q"q") = ( t x )iA tY )ki(t x t Y )^ XY \ (62) 

where X,Y(i,j,k,l,m,n) are the gluon(quark) colour indices and 7 T (XY) represents gluon 
index permutations. The symbols X and Y refers to any combination of internal gluons 
X,Y = 1, ...8 interconnecting the three quark-antiquark lines of colour ij , kl and mn (and 
flavours q, q' and q", respectively). In this case, the calculation of the colour factors is much 
easier. In fact, it is trivial to see that there are only four fundamental ones, of the form 
Ctf, with o,p = 1, 2, 3. Two correspond to the case o = p, 4/3 and —1/6, and two to the 
opposite condition, o/j), —1/3 and 7/6. These values are obtained depending on the way 
three coloured fermion loops (of flavour q ^ q' ^ q") are connected to each other via four 



















Tl ' [Tlijki V '9 9 ) ( T 2 ijki ( qqq'q'gg )) f ] 

/07T,7T / 

°12 

X78 

78X 

8X7 

87X 

X87 

7X8 

X78 

-3(31) 

57/2(-l/2) 

-3(31) 

-3(-5) 

6(-14) 

-3(-5) 

78X 

-3(31) 

-3(31) 

57/2(-l/2) 

-3(-5) 

-3(-5) 

6(-14) 

8X7 

57/2(-l/2) 

-3(31) 

-3(31) 

6(-14) 

-3(-5) 

-3(-5) 

87X 

-3(-5) 

-3(-5) 

6(-14) 

-3(31) 

-3(31) 

57/2(-l/2) 

X87 

6(-14) 

-3(-5) 

-3(-5) 

57/2(-l/2) 

-3(31) 

-3(31) 

7X8 

-3(-5) 

6(-14) 

-3(-5) 

-3(31) 

57/2(-l/2) 

-3(31) 


Table 8 : The colour factors CI 2 , multiplied by 9, for e + e _ —> qqq'q'gg (q 7 ^ q'). In round 
brackets are given the additional ones needed in the case of identical flavours q — q' (for 
which the exchange of indices i ^ k is understood in one of the two interfering colour 
structures), multiplied by 27. Permutations 7r are along columns, permutations 7 r' are along 
rows. The first three are cyclical permutations, whereas the latter three are anticyclical. 
Note the symmetry in the colour factors: the (1 — 3) x (1 — 3) submatrix is identical to the 
(4 — 6 ) x (4 — 6 ) one transposed, whereas the (1 — 3) x (4 — 6 ) submatrix is identical to the 
(4 — 6) x (1 — 3) one transposed. 


Tr [Tiljki ( qqq'q'gg ) {qm'q'gg )) f ] 

^-Y7T,7T / 7r // 

^13 

(X7)(X8) 

(X7)(8X) 

(7X)(X8) 

(7X)(8X) 

X78 

l(-5) 

-7/2(-l/2) 

1(31) 

-7/2(71/2) 

78X 

-7/2(71/2) 

1 (-1/2) 

28(31) 

-8(-5) 

8X7 

28(4) 

-8(-32) 

-7/2(-l/2) 

1(4) 

87X 

-7/2(-l/2) 

1(4) 

-7/2(-5) 

1(40) 

X87 

-8(40) 

28(4) 

l(-5) 

-7/2(-l/2) 

7X8 

1(31) 

-7/2(-5) 

-8(-14) 

28(31) 


Table 9: The colour factors Cyf 77 , multiplied by 9, for e + e _ —> qqq'q'gg (q 7 ^ q'). In 
round brackets are given the additional ones needed in the case of identical flavours q — q' 
(for which the exchange of indices i ^ k is understood in one of the two interfering colour 
structures), multiplied by 27. Permutations 7 r are along columns, permutations n' and n" 
are along rows. 





Tr [U't, (?w V 00 ) (US'" (qw't'gs)) 1 

°33 

(X7)(X8) 

(X7)(8X) 

(7X)(X8) 

(7X)(8X) 

(X7)(X8) 

32 (-32) 

-4(4) 

-4(4) 

l/ 2 (-l/ 2 ) 

(X7)(8X) 

-4(4) 

32(40) 

l/ 2 (-l/ 2 ) 

-4(-5) 

(7X)(X8) 

-4(4) 

l/ 2 (-l/ 2 ) 

32(40) 

-4(-5) 

(7X)(8X) 

l/ 2 (-l/ 2 ) 

-4(-5) 

-4(-5) 

32(31) 


Table 10: The colour factors C 33 ,7r n , multiplied by 9, for e + e _ —> qqq'q'gg (q 7 ^ q'). In 
round brackets are given the additional ones needed in the case of identical flavours q = q' 
(for which the exchange of indices i k is understood in one of the two interfering colour 
structures), multiplied by 27. Permutations 7 r and n' are along columns, permutations tt" 
and n'" are along rows. 


Tr [ T £Jki (qqq'q'gg) ( qqq'q'gg )) f ] 

y^y7T7r ,7T 7T 

( - / 34 

(X8)(X7) 

(X8)(7X) 

(8X)(X7) 

(8X)(7X) 

(X7)(X8) 

-3(31) 

-3(-5) 

-3(-5) 

57/2(-l/2) 

(X7)(8X) 

-3(-5) 

-3(-l/2) 

6(40) 

-3(4) 

(7X)(X8) 

-3(-5) 

6(40) 

-3(-l/2) 

-3(4) 

(7X)(8X) 

57/2(-l/2) 

-3(4) 

-3(4) 

-3(-32) 


Table 11: The colour factors ,n 77 , multiplied by 9, for e + e _ —> qqq'q'gg (q 7 ^ q'). In 
round brackets are given the additional ones needed in the case of identical flavours q — q' 
(for which the exchange of indices i k is understood in one of the two interfering colour 
structures), multiplied by 27. Permutations n and n' are along columns, permutations n" 
and 7 t"' are along rows. 


gluons. If o = p, one of the loops has four attached vectors and the others have two. When 
o 7 ^ p, two loops have three connections and the remaining one only one. (Note that in case 
of the colour diagrams of process ( 0 ) gluon permutations along quark lines are constrained, 
as two (and only two) gluon connections from each of the interfering diagrams must always 
be adjacent in the full colour graph.) Once again, it is rather straightforward to calculate all 
diagrams, and corresponding colour structures, differing only by flavour exchanges, given the 
explicit expressions of the formulae ©"(ED and of the colour matrices in eqs. (j5D|)-(p^), 
in terms of external momenta/helicities and flavours, respectively. 




























Also in case of process (|3j), if two or more flavours in the final state are identical, 
more (spinor and colour) diagrams are needed. For the spinor part, the number of graphs 
doubles (triples) depending on whether two(three) flavours coincide and the additional ones 
(with respect to the case of all different flavours) can again be obtained by simple relabelling 


of the relevant momenta and helicities in the formulae (|49|) - (^ID , as previously explained. 
As for colour diagrams, it can immediately be seen that for the cases q = q' ^ q", q' = q" ^ q 
and q" = q ^ q' one has to compute graphs in which two bubbles are connected by 2, 3 or 
4 gluons with 2, 1 and 0 additional gluons coupled to only one of the loops, respectively. In 
correspondence of the three sets of colour graphs, one gets 3, 2 and 1 independent colour 
factors: i.e., (1/18, 5/9,—4/9), (—7/18,1/9) and (—1/3). Finally, in the case q — q' — q", 
the additional colour factors needed for the calculation of the complete ME for process (||) 
have already been derived, as they appear among the numbers given in Tab. |6]. Indeed, 
one only needs five of those: i.e., 31, —5, —1/2, 71/2 and 4 (all to be further divided by 
27). Also in case of identical flavours the same restrictions as above on the possible gluon 
permutations in the colour diagrams squared apply. 

To appropriately combine the spinor and colour parts of the amplitudes into the matrix 
elements of the form ( |52|) for all processes (|l])-(0) is a tedious but trivial labour that we leave 
as an exercise for the reader, as it would only results here in cumbersome expressions carrying 
no additional information with respect to that already given in the previous formulae and 
tables. This procedure has been implemented in our FORTRAN programs, that we make 
available to the public upon request. 
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e e -> 6 jets at LO at LEP1 



Figure 5: Differential distributions in the following shape variables: thrust (T), oblateness 
(O), C-parameter ( C ), heavy ( Mh ), light {Ml) and difference {Ml,) jet mass for the six- 
(continuous line), five- (dashed line) and four-jet (dotted line) rates at LO, at LEP1. The 
‘star’ points correspond instead to the six-jet rate as obtained by using the HERWIG Monte 
Carlo. The C scheme was used with resolution parameter y = 0.001. Notice that the 
distributions have been normalised to a common factor (one) for readability. 
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Figure 6: Differential distributions in the following shape variables: sphericity ( S ), apla- 
narity (S), planarity (P), minor (m) and major (M) value and D-parameter (D) for the 
six- (continuous line), five- (dashed line) and four-jet (dotted line) rates at LO, at LEP1. 
The ‘star’ points correspond instead to the six-jet rate as obtained by using the HERWIG 
Monte Carlo. The C scheme was used with resolution parameter y = 0.001. Notice that 
the distributions have been normalised to a common factor (one) for readability. 
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Figure 8: The energy spectra Ej of six-jet events at LO in the J scheme for y = 0.001, at 
LEP2 (here, y/s = 175 GeV). Jets are energy-ordered such that E\ > E 2 > ... > E 5 > E 6 : 
j = 1 (continuous line), j = 2 (short-dashed line), j = 3 (dotted line), j = 4 (dot-dashed 
line), j = 5 (fine-dotted line) and j = 6 (long-dashed line). Normalisation is to unity. 
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Figure 9: Subsets of Feynman diagrams contributing to qqq'q'gg production via W + W~. 
For simplicity the leptonic part of the diagrams is not shown. Permutations are not shown 
either. 
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e e -> 6 jets at LO at LEP2 



Figure 10: The distributions in the true invariant mass M % j of events of the type ©-© at 
LO in the C scheme with y = 0.001, at LEP2, for the following combinations of parton pairs 
ij: 12 (continuous lines), 13 (dashed lines), 14 (dotted lines) and 23 (dot-dashed lines). 
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e e -» 6 jets at LO at NLC 
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Figure 11: The total cross section of six-jet events at LO in the C scheme, at NLC, decom¬ 
posed in terms of the three contributions e + e _ —»• qqgggg (continuous line), e + e _ —> qqq'q'gg 
(dashed line) and e + e _ —» qqq'q'q"q" (dotted line). 
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Figure 12: The total cross section of six-jet events at LO in the C scheme, at NLC, de¬ 
composed in terms of the contributions containing 6-quarks. Upper plot: e + e~ —> bbgggg 
(continuous line). Middle plot: e + e~ —> bbbbgg (continuous line), e + e _ —■> bbuugg (dashed 
line) and e + e _ —> bbddgg (dotted line). Lower plot: e + e _ —> 666666 (continuous line), 
e + e _ —> bbbbuu (short-dashed line), e + e~ —> bbbbdd (dotted line), e + e _ —> bbuuuu (dot- 
dashed line) and e + e _ —> bbuudd (long-dashed line); here, the case e + e _ —> bbdddd has not 
been reproduced as it visually coincides with e + e _ —> bbbbdd. In the upper lines mj = 0, in 
the lower lines rrib = 4.25 GeV. 




e e -> 6 jets at LO at NLC 




Figure 13: The distributions in thrust (upper plot) and sphericity (lower plot) for e + e _ —> 
tt —> 6-jet events (continuous line) and for those of the type (d)-© (dashed line) at LO in 
the C scheme with y = 0.001, at NLC. Notice that the distributions have been normalised 
to a common factor (one) for readability. 
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Figure 14: The distributions in the reduced invariant mass rriij (|29|) for events of the type 
(0)-(D at LO in the C scheme with y = 0.001, at NLC, for the following combinations of par- 
ton pairs ij: (12)[23]{35} (continuous lines), (13)[24]{36} (short-dashed lines), (14)[25]{45} 
(dotted lines), (15)[26]{46} (dot-dashed lines) and (16)[34]{56} (long-dashed lines), in the 
(upper) [central] {lower} frame. 











e e 


m jets + n photons at LEP1 



y 


Figure 16: The total cross section for photoproduction at LO in the C scheme, at LEP1, 
for the processes: e + e _ —> qqggg'J and e + e~ —> qqq'q'g^f (m=5, n=l: continuous lines), 
e + e~ —> qqgg 77 and e + e - —> qqq'q' 77 (m=4, n=2: dashed lines), e + e“ —> 9^.9777 ( m= 3, 
n=3: dotted line) and e + e _ —> 997777 (m=2, n=4: dot-dashed line). In the cases m=5, 
11 = 1 and m=4, n=2 the upper lines represent the two-quark contributions whereas the lower 
ones refer to four-quark processes. The following additional cuts have introduced on the 
photons: p? > 5 GeV and | cos# 7 | < 0.85. 





Figure 17: A symbolic representation of the decomposition of a quadruple-gluon vertex into 
three triple-gluon vertices. 



